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ABSTRACT 


This  thesis  investigates  the  feasibility  of  automating 
the  design  of  microprocessor-based  digital  filters.  The 
ability  of  a  prototype  design  system  to  successfully  produce 
filter  realizations  is  tested.  General  filter  structures  and 
programming  algorithms  are  presented.  Shortcomings  in  the 
current  version  of  the  design  system  are  determined. 
Modifications  are  made  as  required  to  support  digital  filter 
realizations.  The  feasibility  of  filter  generation  is 
demonstrated  using  realistic  examples  taken  from  the 
literature . 
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I.  INTRODUCTION 


Advances  in  integrated  circuit  and  microprocessor 
technology  have  allowed  increasingly  diverse  applications  of 
these  devices.  The  decrease  in  their  size  and  power 
requirements  and  the  increased  market  for  them  has  resulted 
in  a  decrease  in  costs  while  yielding  increased  computing 
power.  If  this  trend  continues,  an  even  broader  variety 
applications  can  be  expected.  The  effects  of  increased 
complexity  and  the  decrease  in  component  costs  will  lower 
total  hardware  costs.  However,  as  shown  by  Shooman  [Ref.  1], 
software  costs  can  be  expected  to  remain  high  for  three 
reasons :  new  applications  will  require  new  programs  to  be 
written,  the  replacement  of  older,  existing  computers  with 
newer  versions  will  require  either  completely  new  software 
or,  at  the  very  least,  modifications  to  existing  software, 
and  because  programming  is  highly  labor  intensive,  it  will 
be  strongly  affected  by  inflation.  In  fiscal  year  1980, 
fifty-seven  billion  dollars  were  spent  on  computer  systems. 

Of  this  amount,  thirty-two  billion  dollars,  or  fifty-six 
percent,  was  allocated  to  software  [Ref.  2],  It  can  be 
concluded  that  the  most  costly  expenditures  associated  with 
computer  system  development  and/or  modification  are  software 
related,  and  this  situation  can  be  expected  to  remain 
unchanged.  It  is  therefore  to  our  advantage  to  automate 
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software  production.  The  effort,  in  terms  of  manpower, 
required  to  develop  a  given  software  implementation,  as  well 
as  the  monetary  expense  necessary  to  support  it,  can  then  be 
devoted  to  the  research  and  development  of  new  systems  and 
applications . 

The  application  of  digital  computers  to  the  problems  of 
real  time  control  is  an  increasingly  important  aspect  of 
computer  technology,  and  is  representative  of  the  uses  which 
can  be  expected  with  technological  advances.  The  design  of 
such  systems  has  proven  to  be  a  complex  and  expensive 
undertaking.  The  concepts  which  will  reduce  this  complexity 
and  a  methodology  for  the  automated  production  of  real  time 
control  systems  have  been  developed.  [Ref.  3]  While  no 
computer  exists  which  can  duplicate  the  creative  process  of 
design  illustrated  in  Figure  1,  there  are  elements  of  it 
which  can  be  tasked  to  the  computer  for  accomplishment.  In 
order  to  determine  what  portions  of  the  design  process  can 
be  automated,  the  methodology  used  by  the  designer  must  be 
studied.  The  steps  he  follows  usually  involve  combining 
existing  elements  or  components  in  such  a  way  that  the 
desired  system  is  realized.  The  components  that  are  used  are 
selected  from  a  library  of  such  items  which  is  located 
either  in  the  designer’s  memory  or  is  contained  in  a 
physical  listing.  Because  the  components  contained  in  the 
listing  are  regularly  ordered,  a  computer  can  be  used  to 
maintain  it,  and,  when  given  the  specifications  of  a  desired 


system,  it  can  select  the  necessary  entries.  This  portion  of 
the  design  process  has  been  successfully  automated  [Ref.  4] 
and  is  illustrated  in  Figure  2.  The  proposed  design  tool 
merges  the  ideas  of  automated  software  generation  and 
automated  hardware  production  into  one  system.  The 
description  of  this  system  is  the  subject  of  chapter  two. 

The  original  intent  of  the  design  tool  was  to  automate 
the  production  of  the  software  and  hardware  necessary  for 
the  realization  of  real  time  controllers,  but  if  this  tool 
is  to  be  a  truly  useful  design  aid,  it  must  be  applicable  to 
a  wider  variety  of  problems .  The  implementation  of  digital 
filters  is  suitable  for  a  dedicated  microprocessor  system, 
and  is  typical  of  potential  applications  problems. 

Digital  filtering  has  been  the  subject  of  considerable 
research  during  the  past  fifteen  years.  Various  implemen¬ 
tations  have  been  achieved,  including  hardwired  logic, 
special-purpose  computers,  and  general-purpose  computers. 

With  the  development  of  the  microprocessor,  and  the  recent 
increases  in  wordsize  and  speed,  a  new  alternative  is 
available.  The  application  of  the  proposed  design  system 
to  the  implementation  of  microprocessor-based  digital 
filters  is  the  subject  of  this  thesis. 

The  difference  equation  form  of  the  filtering  function  is 
well-suited  to  microprocessor  realization.  However,  the 
algorithm  used  can  be  expressed  in  a  variety  of  ways. 

Chapter  three  discusses  four  possible  implementations.  The 
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Figure  2.  Automated  Hardware  and  Software  Generation 


mathematical  descriptions  discussed  are  chosen  for  their 
simplicity  and  ease  of  implementation,  but  should  not  be 
construed  as  the  only  available  alternatives.  The 
possibilities  £  'e  virtually  unlimited.  The  filtering 
function  can  be  performed  using  one  of  two  general  methods, 
independent  of  the  algorithm  chosen.  The  first  method 
defines  the  transfer  function  as  the  product  of  the  first  and 
second  order  sections.  This  is  the  familiar  cascade  form  of 
the  digital  filter.  The  alternative  method  is  the  parallel 
form,  which  expresses  the  transfer  function  as  the  sum  of 
first  and  second  order  sections.  [Ref.  5]  The  methods  for  the 
general  implementation  of  each  of  these  applications  within 
the  bounds  of  the  design  system  are  also  discussed  in  chapter 
three.  The  optimization  of  the  solution  is  not  of  concern. 

The  best  realization,  as  determined  by  metrics  such  as  memory 
use  and  chip  count,  is  not  the  current  objective  of  the  design 
system.  If  the  implementation  produced  satisfies  the 
requirements  of  the  problem  in  terms  of  timing  constraints 
and  function,  the  solution  is  acceptable. 

The  application  of  the  design  system  to  problems  other 
than  controller  design  is  important  for  several  reasons.  By 
demonstrating  the  ability  of  the  system  to  generate  success¬ 
ful  realizations  for  a  wide  variety  of  problems,  its  overall 
utility  as  a  design  tool  will  be  proven.  Varying  the  type 
of  implementation  problem  presented  to  the  system  extends 
the  applicability  of  the  design  technique  and  refines  our 


understanding  of  the  system.  By  testing  various  applications, 
inadequacies  in  the  current  implementation  of  the  system  can 
be  detected  and  corrected.  The  attempt  to  generate  digital 
filters  using  the  design  system  revealed  a  number  of  such 
problems  and  the  search  for  these  inconsistencies  became  a 
major  portion  of  the  thesis.  The  problems  found  are  presented 
with  corresponding  solutions  in  chapter  four. 

Chapter  five  describes  the  application  of  the  design 
system  to  realistic  problems.  The  concepts  discussed  in  the 
previous  chapters  are  utilized  to  implement  both  the  cascade 
and  parallel  forms  of  a  fourth  order  digital  filter. 

The  system  described  in  this  thesis  represents  a  useable 
design  tool.  Unfortunately,  the  poding  that  the  user  is 
required  to  produce  in  order  to  achieve  the  desired  hardware 
and  software  realization  is  both  awkward  to  use  and  tedious 
to  generate.  The  addition  of  new  realization  volumes,  as 
well  as  the  maintenance  of  existing  ones,  is  a  difficult 
process  as  well.  The  potential  of  the  design  system  is 
therefore  limited  by  our  ability  to  simplify  the  user  input 
specification  and  provide  a  means  by  which  the  software  and 
hardware  library  may  be  modified  to  accommodate  any  potential 
design  problem.  Recommendations  for  these  and  other 
improvements  are  contained  in  the  sixth  and  final  chapter. 


L6 


II.  SYSTEM  OVERVIEW 


The  current  version  of  the  design  system  has  evolved 
from  the  model  proposed  by  Matelan  [Ref.  6]  and  implemented 
by  Ross  [Ref.  7].  Each  of  these  efforts  used  the  design  of 
real  time  controllers  as  the  problem  model  for  development 
of  the  system.  The  concepts  and  terminology  introduced  by 
Matelan  can  be  found  in  the  current  implementation.  The 
results  of  this  work  are  summarized  here. 

A.  PROBLEM  MODEL 

The  first  aspect  of  the  system  to  be  considered  is  the 
problem  model.  In  order  to  produce  a  successful  realization 
for  any  problem,  it  must  first  be  expressed  in  a  form 
understandable  by  the  design  system. 

1 .  Contingency/Task  Pairs 

Each  problem  can  be  expressed  as  a  collection  of 
contingency/task  pairs,  where  a  contingency  is  defined  as  a 
logical  or  relational  function  of  some  input  variable  or 
variables.  The  associated  task  is  executed  when  the 
contingency  is  satisfied.  Therefore,  the  first  step  in 
developing  a  realization  is  to  express  the  problem  in  terms 
of  contingency/task  pairs. 


High  Level  Problem  DescriDtion 


Implicit  m  the  model  of  the  problem  are  rigid 
constraints  on  the  testing  of  each  contingency  and  the  time 
allowed  for  execution  of  the  corresponding  task.  The 
designer  must  specify  these  real  time  requirements  in  the 
statement  of  the  problem.  Matelan  proposed  a  new  high  level 
design  language  called  a  Control  System  Design  Language,  or 
CSDL,  for  this  purpose.  The  language  consists  of  four 
sections:  identification  section,  environment  section, 

contingency  list,  and  procedures  section.  The  identification 
section  is  a  record  of  the  user  identification  of  the 
problem.  The  environment  section  specifies  the  interface 
between  the  device  and  the  process  which  is  to  be  operated 
upon.  It  defines  all  input  and  output  variables  and  their 
characteristics,  as  well  as  those  variables  internal  to  the 
mathematical  operation  of  the  device.  The  contingency  list 
is  a  declaration  of  those  conditions  that  the  device  must 


respond  to,  the  associated  task  that  each  must  execute,  and 
the  real  time  constraints  imposed  upon  each  pair.  The  timing 
constraints  are  determined  by  the  maximum  time  allowed  to 
recognize  a  contingency  and  the  maximum  time  available  to 
execute  the  corresponding  task.  The  routines  which  comprise 
each  contingency/task  are  contained  in  the  procedures 
section.  By  definition,  contingencies  are  written  as 
functions,  while  tasks  are  specified  as  procedures.  Each 
contains  the  high  level  descriptions  necessary  for  the 


performance  of  its  role  in  the  system  being  produced.  Figure 
3  is  an  example  CSDL  listing  of  a  simple  filtering  problem. 

The  identification  section  is  readily  found  and  easily 
understood.  The  environment  section  specifies  two  variables, 
one  for  input  and  one  for  output.  Each  contains  eight 
TTL-compatible  bits.  The  contingency  section  indicates  that 
the  function  READY  is  executed  every  millisecond,  and  when 
true,  the  task  FILTER  is  performed.  The  procedures  section 
contains  the  descriptions  of  the  steps  necessary  to  perform 
the  test  for  the  contingency  READY  and  the  execution  of  the 
process  FILTER.  In  the  case  of  the  contingency  READY,  an 
external  variable,  DATA,  is  read.  If  its  value  is  one, 
indicating  the  presence  of  data  for  processing,  the 
contingency  is  satisfied  and  the  value  of  READY  is  set  to 
one,  and  the  task  FILTER  is  executed.  If,  on  the  other  hand, 
DATA  is  equal  to  zero,  the  contingency  will  not  be  satisfied 
and  READY  will  remain  equal  to  zero  to  indicate  a  false 
condition.  The  task  FILTER  is  not  executed  under  these 
circumstances . 

The  function  FILTER  contains  the  description  of  the 
difference  equation 

y(n)  =  x(n)  +  0.5y(n-l)  -  0.5y(n-2). 

The  value  of  the  input  is  assumed  to  be  in  digital  form  and 
is  read  first.  The  value  of  the  corresponding  filtered 
output  is  then  computed  and  the  result  is  issued  as  an  output. 


IDENTIFICATION: 

DESIGNER:  M.  R.  Heilstedt 

DATE:  19  Apri 1  1983 

PROJECT:  CSDL  Example  -  Filter 

ENVIRONMENT: 

INPUT:  X,8,TTL 
OUTPUT:  Y , 8/ TTL 

CONTINGENCY  LIST: 

When  READYMMS  Do  FILTER 

PROCEDURES: 

Contingency  READY 
Sense  (DATA)? 

If  DAT  As  1  Then  READY : s 1  fi? 
Exit  READY 

Task  FILTER 

Sense  (X); 

Y=X+(0.5*YN1 )-(0.5*YN2) ? 
Issue  ( Y ) ; 

YN2= YN  1  ; 
ynisy; 

Exit  FILTER 


Figure  3.  Example  CSDL  Listing 


The  current  version  of  the  design  system  does  not 
support  the  use  of  CSDL  for  problem  specifications,  requiring 
that  the  user  generate  the  list  of  primitives  manually.  A 
user  specification  format  based  on  the  Ada  programming 
language  is  currently  under  development.  [Ref.  8]  In  the 
interim,  CSDL  serves  as  a  useful  tool  in  understanding  the 
structure  of  the  problem  specification  required  by  the 
design  system  and  the  transformation  it  undergoes  in  the 
course  of  generating  a  design. 

B.  INTERMEDIATE  FORM  OF  THE  PROBLEM 

The  input  specification  is  translated  into  an  inter¬ 
mediate  representation  consisting  of  a  list  of  primitives. 
This  transformation  is  analogous  to  the  operation  performed 
by  a  compiler  on  high  level  languages  such  as  Fortran  and 
Pascal.  The  list  of  primitives  is  a  complete  description  of 
the  problem,  containing  all  of  the  information  necessary  to 
generate  a  hardware  and  software  implementation  of  the  design 
problem  and  to  verify  that  all  timing  requirements  have  been 
satisfied  with  the  selected  realization  volume. 

A  second  file,  called  the  IADEFL  file,  is  generated  at 
the  same  time  as  the  primitive  list,  and  contains  the  list 
of  contingency/task  pairs  and  their  corresponding  timing 
constraints.  This  data  is  extracted  from  the  ID  table,  the 
environment  table,  the  timing  data  from  the  contingency 
list,  and  the  design  criteria  table.  The  design  criteria 


table  is  generated  from  the  data  listed  in  the  design 
criteria  section  of  the  CSDL  listing.  This  section  was 
added  by  Ross  to  provide  a  method  for  the  designer  to  specify 
the  criteria  upon  which  an  acceptable  realization  could  be 
chosen.  The  criteria  are:  first  realization  generated, 
least  costly  realization,  and  the  realization  that  uses  the 
least  power.  The  current  version  of  the  system  only  supports 
the  first  realization  generated  criterion.  The  creation  of 
the  IADEFL  file  and  the  primitive  list  is  the  first  action 
taken  by  the  design  system  in  its  attempt  to  generate  a 
hardware  realization.  The  operation  of  the  design  system  in 
relation  to  these  files  is  depicted  in  Figure  4. 

C.  FUNCTIONAL  ELEMENTS 

The  problem  formulation  is  performed  by  the  user  for 
both  the  high  level  description  and  intermediate  represen¬ 
tation.  Each  of  these  tasks  will  eventually  be  automated 
as  well,  providing  a  user  friendly  interface  to  the  design 
system.  As  stated  previously,  development  of  the  input 
module  is  in  progress,  but  until  a  satisfactory  input 
specification  to  intermediate  form  translator  can  be 
developed,  the  input  to  the  design  system  must  be  manually 
compiled. 

1 .  Optimizer  Module 

After  the  input  files  have  been  produced,  they  are 
passed  to  the  first  component  in  the  design  system,  the 
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translator 


Figure  4.  Flowchart  of  Automated  Design  System 


optimizer  module.  It  serves  two  purposes,  acting  as  the  main 
program  controlling  the  operation  of  the  system,  and  func¬ 
tioning  as  the  input  module  for  the  primitive  list.  Although 
not  implemented  in  the  current  version  of  the  system,  this 
module  will  have  the  capability  of  comparing  the  multiple 
realizations  generated  by  the  system  and  selecting  that  one 
which  best  satisfies  the  chosen  metric  as  specified  by  the 
user  in  the  design  criteria  section  of  the  input  listing. 

2 .  Functional  Mapper 

The  Functional  Mapper  is  the  first  module  called  by 
the  Optimizer,  and  is  used  to  ensure  that  each  primitive 
contained  in  the  intermediate  problem  specification  can  be 
realized  with  the  current  realization  library.  This  is 
accomplished  using  two  separate  operations.  The  realization 
volume  index  is  first  searched  for  the  primitive  name  as 
given  by  the  current  line  of  the  intermediate  list.  If  the 
primitive  is  found,  the  specifications  associated  with  the 
primitive  are  compared  to  those  contained  in  the  realization 
library.  The  mapping  is  considered  successful  only  if  all  of 
the  primitives  can  be  realized  and  their  criteria  satisfied. 

If  a  given  primitive  cannot  be  mapped  to  a  realization,  or 
if  its  selection  criteria  fail  to  fit  those  of  the  realization 
a  new  realization  library  will  be  searched.  Failure  to 
satisfactorily  realize  a  primitive  in  any  library  results  in 


an  unsuccessful  mapping. 


3 .  Timing  Analyzer 

The  second  module  executed  is  the  Timing  Analyzer. 

It  generates  the  monitor  necessary  for  controlling  the 
operation  of  the  device  under  development ,  ordering  the 
execution  of  the  realization  contingency  tests  and  task 
executions  such  that  the  timing  constraints  will  be  satisfied. 
In  order  to  make  such  an  assurance,  the  timing  analyzer 
assumes  that  all  contingencies  are  true.  This  requires  that 
all  of  the  tasks  must  be  executed  by  the  monitor  and  there¬ 
fore  defines  the  worst  case  situation.  The  resulting  premise 
is  that  if  the  worst  case  can  satisfy  the  timing  requirements, 
all  cases  satisfy  them  as  well.  The  timing  analyzer  determine 
the  length  of  time  required  to  execute  all  of  the  code  con¬ 
tained  in  each  of  the  contingency/task  pairs  and  compares  it 
to  the  timing  constraints  listed  in  the  Applications  Timing 
Table  as  generated  from  the  IADEFL  file.  If  all  of  the 
timing  constraints  are  met,  the  realization  is  successful. 

If  not,  the  realization  fails.  The  output  of  the  Timing 
Analyzer  is  used  to  generate  monitor  primitives  for  success¬ 
ful  realizations.  These  primitives  determine  the  sequence  of 
execution  for  the  contingency/task  pairs,  and  are  added  to 
the  list  of  primitives  derived  from  the  high  level  description 
of  the  problem. 

An  important  result  of  the  research  conducted  by  Ross 
was  that  the  design  system  is  capable  of  automatically 
producing  multi-processor  realizations.  This  is  done  in  the 
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Timing  Analyzer.  If  a  single  processor  realization  is 
impossible,  the  Timing  Analyzer  partitions  the  Applications 
Timing  Table  into  two  separate  lists.  Each  of  these  contains 
the  timing  information  corresponding  to  complete  sets  of 
contingency/task  pairs.  In  other  words,  the  partitioning  of 
the  problem  results  in  groupings  of  these  pairs.  The 
regularly  ordered  nature  of  the  contingency/task  model  makes 
this  possible.  In  order  to  generate  separate  sets  of  hard¬ 
ware,  the  two  timing  tables  are  passed  to  the  Formatter 
module  individually.  Theoretically,  it  is  possible  to  pro¬ 
duce  a  realization  in  which  each  contingency/task  pair  is 
located  in  its  own  processor.  It  is  important  to  note  that 
such  a  system  may  require  shared  resources.  The  current 
version  of  che  design  system  will  only  test  for  a  dual 
processor  realization  when  a  single  one  fails.  The 
generation  of  such  a  realization  cannot  be  forced. 

4 .  Formatter  Module 

The  Formatter  receives  the  complete  list  of 
primitives  necessary  to  produce  the  output  listing  for  the 
design.  The  listing  consists  of  the  software  necessary  to 
perform  the  desired  function  and  the  hardware  required  to 
support  it.  The  hardware  and  software  listings  are  written 
to  separate  files  which  are  in  turn  copied  to  an  output 
device  such  as  a  terminal  or  line  printer.  The  design 
program  is  then  terminated. 


D.  REALIZATION  LIBRARY 

The  realization  library  has  thus  far  been  mentioned  only 
in  describing  the  operation  of  the  other  components  of  the 
design  system,  but  its  importance  in  the  determination  of  a 
successful  design  should  not  be  overlooked.  Each  realization 
library  contains  volumes  of  hardware  and  software  primitives 
based  on  individual  processor  families.  The  only  volume 
currently  implemented  uses  the  Intel  8080  microprocessor  and 
its  various  support  devices  to  generate  realizations. 

The  general  format  of  each  line  in  a  volume  is  the  same. 
Each  line  is  assigned  a  unique  identifying  number,  and 
contains  a  maximum  of  eighty  characters.  The  first  set  of 
lines  within  the  volume  serve  as  an  index  to  the  primitives 
it  contains.  The  lines  of  the  index  are  copies  of  the  title 
line  of  each  entry.  The  current  format  of  the  realization 
volume  allows  a  maximum  of  9999  lines.  The  index  is  not 
considered  when  determining  the  number  of  lines  contained  in 
the  volume . 

Each  line  of  the  volume  must  conform  to  a  specific 
format,  of  whit  there  are  ten:  Primitive  Title  line, 

Comment  line.  Calc  line,  Attr  line,  Call  line,  Include  line, 
If  line,  Begin  Text  line,  End  Text  line,  and  Text  line.  The 
title  line  begins  with  an  S  or  H,  to  denote  hardware  or 
software,  and  contains  the  name  of  the  primitive.  As  used  in 
the  current  version  of  the  design  system,  it  is  the  most 
important  of  the  lines  found  in  the  volume,  providing  the 


correct  format  for  generating  each  of  the  entries  contained 
in  the  intermediate  problem  description.  The  calling 
arguments,  selection  criteria,  and  attributes  of  the 
primitive  are  enclosed  in  parentheses  following  its  name. 

The  attributes,  which  vary  depending  on  the  nature  of  the 
primitive,  define  such  parameters  as  power  consumption, 
latency,  and  chips  used.  Any  or  all  of  the  arguments, 
selection  criteria,  and  attributes  may  be  omitted,  but  the 
commas  that  separate  them  must  appear.  The  comment  line  is 
denoted  by  the  appearance  of  the  letters  C-O-M  as  the  first 
characters  on  the  line.  The  text  that  follows  is  ignored  by 
the  system.  The  Calc  line  allows  the  use  of  global  variables 
within  the  system.  An  example  of  its  application  is  the 
counter  variable  used  to  total  the  number  of  chips  used  in  a 
particular  design.  In  the  current  realization  library,  this 
variable  is  called  ICN  and  is  incremented  by  any  of  the 
hardware  primitives  that  contain  integrated  circuits.  The 
Attr  line  is  similar  to  the  Calc  line,  but  is  used  to 
compute  the  value  of  an  attribute  of  the  primitive  realiza¬ 
tion.  Incl  and  Call  lines  are  used  to  invoke  other  primitives 
from  within  a  primitive.  The  difference  between  the  two  is  in 
the  placement  of  the  output  that  each  generates.  The  output 
from  a  Call  is  inserted  immediately  following  any  previously 
generated  output.  Output  from  an  Incl  statement  will  be  added 
to  that  of  the  primitive  lists  after  all  other  output  from 
the  including  primitive  has  been  produced.  The  If  line 


provides  more  flexibility  in  library  construction  by 
allowing  branch  instructions  within  the  realization  volume. 
The  Begin  and  End  Text  lines  are  used  to  reproduce  descrip¬ 
tive  lines  in  the  files  generated  for  system  output  by  the 
Formatter  module.  These  lines  are  in  the  final  category  of 
text  lines.  The  most  common  use  of  text  lines  is  for 
assembly  code  associated  with  a  software  primitive.  Figure 
5  is  an  example  of  a  short  realization  volume.  The 
construction  of  the  volume  is  further  demonstrated  and 
clarified  in  the  text  of  the  analog  to  digital  converter 
realization  developed  in  chapter  four. 


V2222S. SAMPLE  (P 1 , P2 : 0 , 8, 0 , 5 : 1 0 , 1 0 , - 1 7 , 1 8, 1 9 , 2222 , 225* ) 

V2223C0M  THIS  IS  A  COMMENT  DESCRIBING  S. SAMPLE.  PI  AND 
V2224C0M  P2  ARE  THE  DUMMY  ARGUMENTS  OF  S. SAMPLE.  0  AND  8 
V2225CQM  ARE  THE  MINIMUM  AND  MAXIMUM  VALUES  OF  THE  ACTUAL 
V2226C0M  ARGUMENT  REPRESENTED  BY  PI,  0  AND  5  ARE  THE 
V2227C0M  CORRESPONDING  MINIMUM  AND  MAXIMUM  FOR  THE  ACTUAL 
V2228C0M  ARGUMENT  REPRESENTED  BY  P 2.  THE  10,10  INDICATES 
V2229C0M  THAT  THE  STORAGE  AND  TIME  ATTRIBUTES  FOR  THIS  ENTRY 
V2230CQM  ARE  EACH  OF  VALUE  10.  THE  -17  INDICATES  THAT  THE 
V2231C0M  VALUE  OF  THE  EXTERNAL  ATTRIBUTE  IS  CALCULATED  ON 
V2232C0M  LINE  2239  C 2222- ( - 1 7) =2239)  .  THE  18  INDICATES  THAT 
V2233COM  LINE  22«0  (2222+ 1 8=22«0 )  CALLS  FOR  CALCULATION  OF 
V2234C0M  SOME  GLOBAL  VALUE.  THE  19  INDICATES  THAT  LINE  2241 
V223SCOM  CALLS  FOR  THE  INCLUSION  OF  SOME  OTHER  REALIZATION. 
V2236C0M  2222  AND  2258  ARE  THE  FIRST  AND  LAST  LINE  NUMBERS 
V2237COM  OF  THIS  REALIZATION. 

V2238C0M 

V2239ATTR  EXTERNAL  =  DUMMY  1  *  7 
V2240CALC  TOTEVT  =  TOTEVT+1 
V2241INCL  H . ALSO SAMPLE (0,0) 

V2242C ALL  S . SAMPTHREE  (TOTEVT) 

V2243C0M 

V2244C0M  H.ALSOSAMPLE  AND  S. SAMPTHREE  ARE  TwO  OTHER  ENTRIES. 
V2245C0M  ALSOSAMPLE  HAS  TWO  DUMMY  ARGUMENTS,  BOTH  ASSIGNED 
V2246C0M  VALUE  ZERO  IN  THIS  CASE.  S. SAMPTHREE  HAS  ONE  DUMMY, 
V2247C0M  SET  EQUAL  TO  THE  GLOBAL  TOTEVT  FOR  THIS  CASE. 
V2248C0M 

V2249IF  DUMMYl  .GE.  2  SKIP  2 
V2250INCL  H.SAMPFOUR  (  ) 

V2251C0M  INCLUDE  H.SAMPFOUR  ONLY  IF  DUMMY!  IS  LESS  THAN  TWO. 
V2252C0M  BEGIN  THE  TEXT  BLOCK  AFTER  THE  NEXT  LINE.  . 
V2253BEGIN  STEXT 

V2254  EVERYTHING  IN  THIS  BLOCK  BETWEEN  BEGIN  AND  END  IS 

V2255  COPIED  TO  THE  OUTPUT  LISTING  EXCEPT  FOR  DUMMY 

V2256  ARGUMENTS  AND  GLOBALS  WHICH  ARE  REPLACED  BY  THEIR 

V2257  ACTUAL  ARGUMENTS  OR  VALUES. 

V2258ENOTEXT 


Figure  5.  Sample  Realization  Library  Entry 
[Ref.  8] 
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III.  FILTERING  METHODS 


With  the  advent  of  the  microprocessor,  a  new  and 
interesting  device  for  the  implementation  of  digital  filters 
has  been  created.  Increases  in  word  length  and  computational 
speed  will  encourage  more  complex  filtering  applications 
using  this  device.  It  is  therefore  important  to  investigate 
the  feasibility  of  producing  digital  filter  realizations 
using  the  automated  design  system.  The  fundamental  require¬ 
ments  that  the  application  places  on  the  system  can  be 
determined  using  the  Intel  8080  microprocessor  library 
currently  available.  Despite  the  filtering  limitations 
presented  by  its  eight  bit  format,  the  deficiencies  present 
in  the  current  system  can  be  identified  and  corrected.  The 
resulting  improvements  will  provide  a  more  useful  design 
tool  and  make  possible  the  generation  of  satisfactory 
solutions  to  more  demanding  problems  when  additional 
realization  volumes  constructed  around  advanced 
microprocessors  are  added  to  the  library. 

A.  TRANSFER  FUNCTION 

The  general  form  of  the  digital  filter  transfer  function 


where  a^  and  b^  are  real  coefficients  and  n  is  the  maximum  of 
the  orders  of  the  numerator  and  denominator  polynomials . 

Figure  6  is  a  flow  diagram  representation  of  this  transfer 
function.  There  are  a  number  of  similar  structures  that  can 
be  used  to  describe  the  basic  filtering  operation.  Those  in 
which  the  real  coefficients  a^  and  b^  are  multipliers  are 
referred  to  as  direct  structures.  Nagle  and  Nelson  [Ref.  10] 
describe  four  direct  structures  and  the  equations  associated 
with  them. 

1 .  First  Direct  Structure 

The  first  direct  structure  is  derived  from  the 

equation 

n  n  . 

H(z)  :  (  I  a.  z~x)/(  E  b.  z~x). 
i=0  1  i=0  1 

The  coefficient  bQ  is  equal  to  one.  Defining  the  input  to  the 
filter  as  X(z)  and  the  output  as  Y(z),  the  previous  equation 
can  be  rewritten  as 


Y(z) 

xTzT 


n  n 

(  E  a.  z’x)/(  e  b.  z~x). 
i=0  1  i=  0  1 


This  can  be  redefined  using  an  intermediate  variable  M(z), 
such  that 


Y(z)  M(z) 
RT7T  *  Ym 


n 

(  E  a. 
1 


z"1)/ 


n 


b. 


z_x). 


[( 


The  following  relationships  can  now  be  derived: 


Y(z) 

MTzT 


n 


) 


and 


M(z)  _  1 

X( z)  “  n 

2  b.  z"1 
i  =  0  1 


or 


X(z) 

MTzT 


n 
:  2 
i=  0 


b. 

i 


). 


Rearranging  these  equations  provides  expressions  for  the 
input  X(z)  and  the  output  Y(z)  in  terms  of  M(z): 


n 

Y(z)  s  (  2  a.  z"1)H(z) 
i  =  0  1 


and 

n 

X(z)  =  (  2  b.  z_1)M( z) 


Because  the  input  is  a  known  quantity,  it  is  desirable  to 
derive  an  equation  for  M(z)  in  terms  of  X(z).  Manipulation 
of  the  previous  equation  for  X(z)  gives 

n 

M(z)  =  X(z)  E  b.  z-1M(z) 
i=l  1 

Because  the  calculations  necessary  to  produce  the  filtering 
function  are  to  be  performed  in  real  time,  the  inverse 
z-transform  of  M(z)  and  Y(z)  is  taken  to  provide  a  time 
domain  solution  of  the  first  direct  structure.  The  results 
are: 


n 

m(k)  =  x(k)  E  b.  m(k-i) 

i  =  l  1 


and 


n 

y(k)  =  E  a.  m(k-i). 
i  =  0  1 

The  flow  diagram  of  the  first  direct  structure,  of  order 
two,  is  shown  in  Figure  7. 

2.  Second  Direct  Structure 

The  transpose  of  a  digital  filter  structure  is 
formed  by  reversing  the  signal  flow  in  all  of  the  branches 
of  the  block  diagram.  Its  transfer  function  is  the  same  as 
that  of  the  structure  from  which  it  was  derived.  The  second 
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direct  structure  is  generated  by  applying  this  principle  to 
the  first  direct  structure.  The  corresponding  time  domain 
solution  also  implements  the  same  filtering  function,  but  an 
additional  difference  equation  is  required.  The  equations 
for  the  second  direct  structure  are: 


p^(k)  =  p^+1(k-l)  +  a^x(k)  -  b^(y(k));  i=l,...,n-l 

p  (k)  =  a  x(k)  -  b  y(k) 
n  n  nJ 

y(k)  =  a^x(k)  +  p^(k-l). 


The  second  direct  structure,  of  order  two,  is  illustrated  in 


the  flow  diagram  of  Figure  8 . 
3.  Third  Direct  Structure 


To  obtain  the  third  direct  structure,  the  expression 


H(z) 


Y(z) 

xTzT 


n  .  n 

(  Z  a.  z~x)A  Z 
i  =  0  1  i=  0 


b. 

l 


is  rewritten  as 


n 

Y(z) (  Z  b.  z-1) 
i  =  0  1 


n 

X(z)(  Z 
i=0 


a . 
i 


z”1) . 


The  output  expression  is  then 


Y(z)  = 


n 

z 

i  =  0 


z-1X( z)  - 


n 

Z 

i  =  l 


b. 

i 


z_1Y(z) . 
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Once  again  taking  the  inverse  z-transform,  the  time  domain 
solution  is  derived: 

n  n 

y(n)  =  Z  a.  x(n-i)  -Zb.  y(n-i). 
i=0  1  i=l  1 

Figure  9  illustrates  the  flow  diagram  for  the  third  direct 
structure  of  order  two. 

4 .  Fourth  Direct  Structure 

The  fourth  direct  structure  is  derived  by  taking  the 
transpose  of  the  third  direct  structure.  The  time  domain 
representation  is: 

rQ(k)  =  x(k)  +  r1(k-l) 

ql(k)  =  anr0(k) 

r  (k)  =  -b  r  (k) 
n  n  0 

q^(k)  =  a^rQ(k)  +  q^+^(k-l),  i=l,...,n-l 
r.(k)  =  -b.rn(k)  +  r.,,(k-l). 

l  l  0  i+l 

The  flow  diagram  for  the  fourth  direct  structure,  of  order 
two,  is  shown  in  Figure  10. 

B.  SECOND  ORDER  SECTIONS 

Each  of  the  structures  described  applies  to  the  general 
case  in  which  the  filter  transfer  function  is  of  order  n. 

The  filtering  function  can  be  successfully  realized  using 
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any  of  these  methods,  provided  that  n  does  not  exceed  two. 

For  transfer  functions  of  higher  order,  coefficient 
sensitivity  becomes  a  significant  factor  in  accurately 
generating  the  filter  output.  This  phenomenon  is  created  by 
the  effects  of  quantization  errors  in  the  representation  of 
the  transfer  function  coefficients.  The  filter  structures 
and  their  related  difference  equations  assume  finite 
precision  in  the  filter  parameters.  However,  because  the 
word  length  used  by  microprocessors  is  fixed,  the  precision 
with  which  the  filters  can  be  realized  is  limited.  This 
means  that  the  coefficient  representations  are  not  exact, 
resulting  in  a  set  of  poles  and  zeros  for  the  realization 
that  differ  from  those  desired.  The  frequency  response  of 
the  actual  filter  will  therefore  be  different  from  the 
design  specification.  In  the  most  extreme  case,  the  filter 
will  become  unstable  if  one  or  more  of  the  poles  are  moved 
outside  of  the  unit  circle  [Ref.  11],  As  shown  by  Rader  and 
Gold  [Ref.  12],  high  order  filters  can  be  redefined  as 
combinations  of  first  and  second  order  sections  in  order  to 
minimize  these  effects. 

C.  COMBINATIONAL  STRUCTURES 
1.  Cascade  Structure 

The  first  combinational  structure  to  be  considered 
is  the  cascade  form  of  the  digital  filter.  It  is  illustrated 
in  Figure  11(a)  and  is  obtained  by  factoring  the  transfer 


function  into  a  product  of  second  order  sections.  This  is 
represented  by  [Ref.  13] 


M 

H(z)  =  n  H. (z) 
i=l  1 


where 


H.  (z) 
1 


1  +  b 


li 


+  a2i 


+  b 


2i 


For  the  fourth  order  transfer  function 


H(z) 


-1  -2  -3  -2 

aQ  +  a^  z  +  a2  z  +  a3  z  +  a^  z 


- n - r? - n - 

1  +  b,  z  +  b_  z  +  b„  z  +  b  z 

X  i  J  T 


—4 


an  equivalent  representation  is 


H(z)  = 


a01  +  all  Z  1  +  a12  Z  a02  +  al2  2  +  a21  Z 


1  *  bll  2-1  +  b21 


1  *  b12  z'1  ♦  b22  2 


-2 


The  first  element  in  the  cascade  receives  its  input  from  the 
signal  to  be  filtered,  while  each  of  the  remaining  elements 
acts  upon  the  output  of  its  predecessor. 

In  z-transform  notation,  the  cascade  structure  output 


is  written 


where 


Yx(z)  =  H1(z)X(z) 

Y.(z)  =  H.(z)Y.  i ( z )  for  i=2 ,3 , . . . ,M-1 . 

The  equations  necessary  for  a  real  time  implementation  are 
formed  by  taking  the  inverse  z-transform  of  these  expressions. 
For  a  system  of  second  order  modules,  the  general  equation 
for  each  section  is 

y(n>  a0M*M-l(n)  *  +  a2M^M-l(n-^) 

-  b1My(n-l)  -  b2My(n-2) 

where 


The  subscript  i  corresponds  to  the  number  of  modules  needed 
to  realize  the  desired  filtering  function,  and  is  equal  to 
|*n/2"J,  where  n  is  the  order  of  the  denominator  polynomial. 
For  the  case  where  n  is  odd,  one  of  the  modules  will  be  of 
order  one. 
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2.  Parallel  Structure 


The  parallel  form  of  the  digital  filter  is  the 
second  structure  to  be  considered,  and  is  illustrated  in 
Figure  11(b).  It  is  formed  by  factoring  the  transfer 
function  such  that  [Ref.  14] 


H(z)  =  E  H. (z) 
i=l  1 


where 


H(z)  = 


an .  +  a, .  z 
_ Oi  li _ 

1  +  bn  .  z“  +  b0 .  z" 
li  2i 


For  example,  the  fourth  order  transfer  function 


can  be  redefined  as 


Each  second  order  section  receives  the  same  input.  The 
output  of  the  system  is  found  by  summing  the  outputs  of  each 
of  the  parallel  elements  in  the  realization,  and  is 
represented  by 


6 


Y(z) 


M 


The  corresponding  time  domain  equation  is 

y(n)  = 

yl(n)+y2 

(n)+. . ,+y^(n) 

where 

Y^n) 

=  aQ1x(n) 

+  a^x(n-l)  -  b^y(n-l) 

-  b2ly(n-2) 

y2(n) 

=  aQ2x(n) 

+  a12x(n-l)  -  b21y(n-l) 

-  b22y(n-2) 

y^n) 

=  a0.x(n) 

+  a^^xCn-l)  -  b^yCn-l) 

-  b^2y(n-2). 

The  equations  presented  for  the  parallel  and  cascade 
structures  are  sufficient  to  develop  the  high  level  problem 
description  necessary  to  realize  microprocessor-based 
digital  filters  using  the  design  system  described  in  the 
previous  chapter. 

D.  DESIGN  CONSIDERATIONS 

The  current  implementation  of  the  design  system 
generates  the  control  code  governing  execution  of  the 
realization  using  a  polled  monitor  strategy.  Pollack 
[Ref.  15]  argued  that  this  approach  was  not  adequate  for 
general  controller  applications  and  proposed  the  introduction 
of  an  interrupt  driven  monitor  as  a  user  selectable 


alternative.  The  existing  realization  volume  has  the  ability 
to  recognize  an  interrupt,  but  its  immediate  response  is 
limited  to  the  setting  of  a  corresponding  condition  flag. 

The  status  of  the  flag  is  examined  and  the  appropriate  task 
executed  as  required  on  a  time  schedule  basis  which  conforms 
to  the  structure  of  the  other  tasks.  The  theoretical 
development  necessary  for  implementing  an  interrupt  driven 
monitor  has  been  completed  [Ref.  16]  and  will  make  possible 
the  realization  of  a  wider  variety  of  control  devices. 

The  polled  monitor  is  the  preferred  strategy  for  digital 
filter  implementations.  Figure  12  is  a  flowchart  represen¬ 
tation  of  three  operations  basic  to  the  filtering  function. 
Each  of  these  tasks  is  performed  once  for  each  sample  taken. 
The  interval  of  time  which  elapses  between  successive 
executions  of  each  contingency/task  pair  is  determined  by 
the  sampling  interval  required  for  the  signal  being  pro¬ 
cessed.  Only  one  execution  is  allowed  during  that  time. 
Therefore,  the  generation  of  the  filter  function  is  periodic, 
with  one  output  calculated  each  period  of  the  sampling 
frequency.  Because  the  polled  monitor  results  in  the 
execution  of  the  contingency/task  pairs  at  regularly 
determined  intervals,  it  is  well  suited  for  applications  in 
which  periodicity  is  required. 

1.  Sampling  Frequency 

The  Nyquist  criterion  determines  the  rate  at  which 


the  signal  to  be  processed  must  be  sampled  in  order  to 


accurately  produce  the  corresponding  filter  output.  The 
amount  of  time  available  to  perform  the  required  computa¬ 
tions  is  then  equal  to  the  elapsed  time  between  successive 
samples.  This,  in  turn,  is  equal  to  the  period  of  the 
sampling  frequency.  This  means  that  the  total  time  required 
to  test  for  all  contingencies  and  execute  their  corresponding 
tasks  must  be  less  than  or  equal  to  the  sampling  period  in 
order  to  produce  a  successful  single  processor  realization. 

If  this  requirement  cannot  be  met,  the  design  system  will 
partition  the  problem  into  two  sets  of  contingency/task 
pairs  and  attempt  to  generate  a  solution  using  two  pro¬ 
cessors.  In  this  case,  the  contingency/ task  pairs  assigned 
to  a  given  processing  element  must  be  able  to  complete 
execution  during  one  period  of  the  sampling  frequency.  The 
maximum  frequency  that  can  be  filtered  by  a  given  realization 
is  therefore  limited  by  the  speed  with  which  the  filtering 
algorithm  can  be  executed,  which  in  turn  is  a  function  of  the 
microprocessor  family  used  to  generate  the  realization. 

2 .  Structure  Selection 

The  parallel  and  cascade  structures  are  both  readily 
adapted  to  microprocessor  implementation.  Because  both 
methods  are  combinations  of  first  and  second  order  sections, 
it  is  a  straightforward  procedure  to  implement  each  section 
as  a  cont ingency/task  pair.  The  natural  partitioning  that 
exists  in  each  of  these  structures  is  well-suited  to 
implementation  using  the  automated  design  system  and  is 


highly  compatible  with  the  partitioning  performed  for 
multi-processor  realizations. 

The  selection  criterion  to  be  considered  when  choosing 
a  filtering  structure  concerns  the  delay  time  incurred  in 
producing  an  output  from  a  specific  input.  The  delay  time 
will  be  a  constant  value  for  the  parallel  realization, 
independent  of  the  number  of  second  order  modules  required 
to  perform  the  filtering  operation.  This  will  correspond  to 
an  interval  equal  to  twice  the  period  of  the  sampling 
frequency,  assuming  a  multi-processor  realization.  The  delay 
time  associated  with  the  cascade  implementation  is  dependent 
upon  the  number  of  second  and  first  order  modules  contained 
in  the  realization.  The  delay  time  will  then  be  equal  to  the 
number  of  modules  in  the  cascade  times  the  period  of  the 

sampling  frequency.  Therefore,  as  the  order  of  the  filter  ] 

increases,  the  number  of  second  order  modules  required  to 
implement  it  will  increase,  causing  a  corresponding  rise  in 
the  filter  delay  time.  These  relationships  are  illustrated 
in  Figure  13. 

The  effects  related  to  ordering  of  the  second  order 
sections  within  a  cascade  structure  are  not  of  concern  for 
the  application  presented  in  this  thesis,  but  in  practice, 
consideration  must  be  given  to  the  changes  that  can  occur  in 
the  limit  cycle  and  quantization  noise  properties  of  the 
filter  as  the  second  order  modules  are  rearranged  in  the 
cascade.  [Ref.  17] 

> 
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IV. 


MODIFICATIONS  TO  THE  SYSTEM 


Applications  problems  are  important  to  the  development 
of  the  design  system  because  each  one  possesses  specific 
requirements  for  implementation.  The  testing  of  widely 
varying  realizations  will  indicate  those  areas  in  which  the 
system  is  deficient  and  identify  hardware  and  software 
primitives  that  are  required  but  abset  from  the  available 
library.  The  digital  filtering  problem  revealed  several 
shortcomings  in  the  current  implementation  of  the  design 
system. 

A.  PARALLEL  PROCESSING 

1 .  Single  Contingency/Multiple  Tasks 

The  parallel  filter  structure  presented  in  chapter 
three  consists  of  a  finite  number  of  elements  that  act 
concurrently  on  the  same  data.  The  output  of  each  processor 
is  then  passed  to  a  summing  element  which  generates  the 
filter  output.  In  terms  of  the  original  problem  model,  this 
corresponds  to  a  single  contingency  (the  availability  of 
data  for  processing)  and  multiple  related  tasks  (each  of  the 
parallel  processing  elements  and  the  summing  junction) .  The 
design  system  implementation,  as  well  as  the  specification 
language  CSDL,  provide  no  means  for  expressing  such  a 
problem  directly. 


2.  FORK  Construct 


Representation  of  the  single  contingency/multiple 
task  structure  subject  to  the  limitations  of  the  design 
language  can  be  solved  by  creating  a  new  construct.  This 
statement  has  been  added  to  the  syntax  of  C SDL  and  is  called 
FORK.  FORK  allows  the  concurrent  execution  of  two  tasks 
followed  by  a  third,  which  in  turn  requires  data  generated 
by  the  first  two.  The  structure  has  the  general  form 

If  MAIN  then  begin 

Do  FORK (T SKI ,TSK2 ) ; 

Do  TSK3 ; 

End 

which  is  stated  in  the  syntax  of  CSDL  as 

When  MAIN:  TIME  Do  F0RKCTSK1 ,TSK2 ,TSK3) . 

TSK1  and  TSK2  are  the  procedures  to  be  executed  in  parallel 
and  are  referred  to  as  the  "forked"  tasks.  Because  it 
combines  the  results  of  TSK1  and  TSK2  in  a  predetermined 
fashion,  TSK3  is  termed  and  "joined"  task.  The  CSDL 
expression  of  the  FORK  construct  is  not  sufficient  to  permit 
the  use  of  the  single  contingency/multiple  task  structure  in 
design  specifications  because  no  translator  exists  to 
produce  the  intermediate  form  of  the  instruction.  However, 
it  is  possible  to  manually  generate  the  required  high  level 
description  of  the  FORK  construct  and  from  it  determine  the 
intermediate  representation.  To  do  so  requires  that  a  set  of 
"dummy"  contingencies  corresponding  to  each  of  the  tasks  be 
generated.  These  contingencies  test  unique  flags  that  are 
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set  by  an  additional  contingency /task  pair  which  controls 
the  execution  of  the  composite  FORK  structure.  Figure  14  is 
a  general  CSDL  representation  of  the  FORK  construct. 


3 .  Timing  Requirements 

The  FORK  construct  makes  possible  the  implementation 
of  design  specifications  incorporating  single  contingency/ 
multiple  task  algorithms.  The  determination  of  the  timing 
constraints  associated  with  the  "dummy”  contingencies  it 
generates  presents  a  new  problem.  There  are  two  potential 
solutions  for  consideration.  Each  of  these  can  be  related 
to  the  single  and  multiple  processor  realizations  that  the 
design  system  is  capable  of  generating.  The  former  is  the 
more  straightforward  derivation  and  will  be  presented  first. 
The  contingency  and  task  names  of  Figure  14  are  used  to 
provide  clarification. 

Implicit  in  both  developments  is  the  assumption  that 
the  user  will  only  specify  the  time  constraint  associated 
with  the  test  for  the  contingency  MAIN.  The  timing  require¬ 
ments  for  the  dummy  contingencies  Cl,  C2,  and  C3  must  be 
determined  from  this  specification.  The  scheduling  of  the 
test  for  MAIN  determines  how  often  the  concurrent  processes 
must  be  tested.  Therefore,  the  scheduling  interval  of  this 
contingency  will  also  be  used  for  Cl  and  C2.  This  value  will 
be  denoted  as  p .  The  scheduling  requirement  for  contingency 
C3  is  related  to  the  availability  of  data  from  tasks  TSK1 
and  TSK2.  The  worst  case  condition,  that  is,  the  minimum 


Contingency  List? 


When 

MAIN 

:  1 1 

SE 

Do 

FORK 

When 

Cl 

:t2 

SE 

Do 

TSK  1 

When 

C2 

:t2 

SE 

Do 

TSK2 

When 

C  3 

:  t  3 

SE 

Do 

TSK  3 

Procedures? 

Function  MA I N 

8i nary  MAIN,  1  ? 

If  G0=1  Then  MAIN:=1? 

Exit  MAIN 

Function  Cl 

Binary  C  1,1? 

If  VARl=l  Then  Cl :=1 ? 

VAR1 :=0? 

Exit  Cl 

Function  C2 

Bi nary  C2,  1  ? 

If  VAR2=0  then  C2:=l? 

VAR2:=0? 

Exit  C2 

Function  C3 

Bi nary  C3,  1  ? 

If  DONE  1  =  1  .and.  D0NE2=1  then  C3:  =  l? 
DONE  1 J  =0  ? 

D0NE2 : =0  ? 

Exit  C3 


Figure  14.  CSDL  Listing  of  FORK  Construct 
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Task  FORK 

VARi:=l? 

VAR2:=l; 

Exit  FORK 

Task  TSK1 

(generate  outDut#  OUT1) 
OONE1 :  =  1  ; 

Exit  TSK 1 

Task  TSK2 

(generate  outout»OUT2) 
DONE 2: =2? 

Exit  TSK2 

Task  TSK3 

0UT3:=0UT1+0UT2? 

Exit  TSK3 


Figure  14.  (continued) 
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interval  at  which  C3  must  be  tested,  exists  when  Cl  and  C2 
are  always  true.  If  these  contingencies  are  scheduled  for 
testing  every  p  seconds,  new  data  will  be  generated  by  TSK1 
and  TSK2  with  the  same  frequency.  This  implies  that  TSK3 
must  be  executed  at  intervals  of  p  seconds  as  well  in  order 
to  maintain  proper  sequencing  of  the  device.  Therefore,  the 
timing  constraint  specified  for  contingency  MAIN  will  also 
be  applicable  to  Cl,  C2,  and  C3.  An  important  fact  that  must 
be  considered  in  reaching  this  conclusion  is  that  contin¬ 
gencies  Cl  and  C2  can  only  occur  as  often  as  contingency 
MAIN.  Similarly,  C3  can  only  be  true  when  Cl  and  C2  have 
been  true  as  well.  The  timing  diagram  associated  with  the 
scheduling  of  the  FORK  construct  for  the  single  processor 
realization  is  illustrated  in  Figure  15. 

The  scheduling  for  the  multiprocessor  realization  is 
not  as  readily  determined.  As  in  the  serial  case,  contingen¬ 
cies  Cl  and  C2  must  be  tested  at  intervals  of  p  seconds,  the 
timing  specification  of  the  MAIN/FORK  contingency /task  pair. 
The  scheduling  of  C3  is  made  difficult  by  the  possibility 
that  each  of  the  contingency/task  pairs  may  reside  in 
different  processors  and  no  synchronization  exists  between 
these  elements.  If  C3  is  tested  with  the  same  frequency  as 
Cl  and  C2,  it  is  possible  to  ’’lose"  a  set  of  data  from  TSK1 
and  TSK2 .  This  occurs  when  these  tasks  generate  data  for 
use  by  TSK3  near  the  beginning  of  their  respective  executions 
and  TSK3  processes  the  data  near  the  end  of  its  computations. 
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Figure  15.  Timing  Diagram  of  FORK  Construct,  Single  Processor 


If  TSK1  and  TSK2  begin  another  processing  cycle  and  produce 
new  data  before  TSK3  has  completed  its  calculations  on  the 
previous  set  of  data,  its  output  will  be  determined  using 
the  second  data  set  generated  and  the  first  values  will  never 
be  processed.  The  timing  diagram  for  such  an  occurrence  is 
shown  in  Figure  16 . 

There  are  two  possible  solutions  to  prevent  the  loss 
of  a  data  set.  The  first  is  to  decrease  the  interval  of 
testing  for  contingency  C3.  Checking  for  the  presence  of  data 
from  TSK1  and  TSK2  with  the  correct  frequency  will  ensure  that 
it  is  processed  by  TSK3  almost  immediately  after  it  is 
generated.  This  solution  results  in  a  new  question  to  be 
answered:  how  must  of  a  decrease  in  the  scheduling  interval 

is  sufficient  to  prevent  data  loss?  The  answer  is  not 
general  in  nature  because  the  timing  specification  required 
will  be  dependent  upon  the  particular  application  being 
considered.  The  incorporation  of  such  a  construct  in  the 
design  system  would  require  a  more  detailed  analysis  by  the 
design  engineer.  This  is  contrary  to  the  original  goals  of 
the  system  and  is  therefore  unacceptable.  An  alternative 
approach  which  preserves  application  independence  is  to  mark 
the  data  as  it  is  produced  by  TSK1  and  TSK2  and  save  it  for 
further  processing  by  TSK3 .  This  implies  the  creation  of  a 
queue  in  which  the  data  is  placed  to  await  action  by  TSK3. 

The  implementation  of  this  proposal  requires  two  queues,  one 
each  for  TSK1  and  TSK2.  Each  time  a  task  generates  new  data, 


Figure  16.  Data  Loss  in  Multiprocessor  FORK  Construct  Implementation 


a  counter  associated  with  its  queue  is  incremented  to  record 
the  entry  of  a  new  element  in  the  service  line.  The 
contingency  C3  then  tests  the  value  of  each  counter.  If  both 
are  greater  than  zero,  the  data  that  has  been  in  each  queue 
the  longest  is  removed  for  processing.  The  process  is 
illustrated  in  Figure  17. 

The  creation  and  maintenance  of  these  data  queues 
solves  the  problem  of  scheduling  C3.  The  contingency  need 
only  be  tested  at  intervals  corresponding  to  p.  If  the  queue 
counters  are  greater  than  zero,  data  has  been  produced  and 
is  ready  for  processing.  The  criteria  upon  which  the 
existence  of  a  true  condition  for  contingency  C3  is  based 
has  therefore  been  changed  from  the  generation  of  output  by 
TSK1  and  TSK2  to  the  presence  of  data  for  processing  in 
their  respective  queues.  This  approach  to  scheduling  of 
the  FORK  construct  does  not  necessitate  the  maintenance  of 
separate  algorithms  for  the  serial  and  parallel  case.  The 
use  of  data  queues  to  insure  accurate  output  from  a  device 
is  suitable  for  single  processor  realizations  as  well. 

Because  the  execution  of  contingency/task  pairs  is  synchronized 
for  such  designs,  data  is  processed  by  TSK3  following  its 
generation  by  TSK1  and  TSK2 .  The  output  of  each  of  these 
tasks  is  stored  in  its  respective  queue  and  is  then  removed 
for  processing  by  TSK3.  Therefore,  the  number  of  individual 
entries  in  either  queue  will  not  exceed  two.  The  CSDL 
specification  of  the  modified  FORD  construct  is  shown  in 
Figure  18. 


One 


g  Diagram  of  FORK  Construct,  Multiple  Processors 


Contingency  List* 

When  MAINstl  SE  Do  FORK 
When  Cl  :t2  SE  Do  TSK1 
When  C2  :t2  SE  Do  TSK2 
When  C3  : t  3  SE  Do  TSK  3 

Procedures  * 

Contingency  MAIN 

Bi nary  MAIN*  1  * 

If  G0=1  Then  MAIN:=l; 

Exit  MAIN 

Contingency  Cl 

Binary  C 1  *  1  * 

If  V AR 1  - 1  Then  Cl:  =  t,* 

VAR  1 S  =0 ; 

Exit  Cl 

Contingency  C2 

Binary  C2»  l  J 
If  VAR2=0  then  C2:=i; 

VAR2:=o; 

Exit  C2 

Contingency  C3 

Binary  C3* 1  * 

If  CTRl.gt.O  .and.  CTR2.qt.O  then  C3:=l* 
Exit  C3 


Figure  18.  CSDL  Listing,  Modified  FORK  Construct 


Task  FORK 

VAR1 :=1 ? 

VAR2:=1? 

Exit  FORK 

Task  TSK 1 

(generate  output#  OUT1) 
Q( 12) :=OUTl ; 

CTR1 : =CTR1 ♦ 1 # 

Exit  TSK  1 

Task  TSK2 

(generate  outout#0UT2) 
Q(22):=OUT2? 

CTR2:=CTR2  +  1  #* 

Exit  TSK2 

Task  TSK3 

Q(tl):=QU2); 

0(21) :sQ(22) ; 

0I.IT3:=Q(  l  1  )+Q(2l )  ; 

CTR1 jsCTRl-1 ; 
CTR2:=CTR2-1  #* 

Exit  TSK 3 


Figure  18.  (continued) 


MULTIPLE  REFERENCES  TO  I/O  VARIABLES 


1 .  Hardware  Binding 

The  organization  of  the  Intel  8080  realization 
volume  as  implemented  in  the  current  design  system  is  based 
on  the  concept  of  hardware  binding  developed  by  Matelan 
[Ref.  13].  It  specifies  that  all  of  the  software  required 
for  a  particular  implementation  must  be  generated  first, 
followed  by  the  hardware  necessary  to  support  it.  The 
present  version  of  the  design  program  modifies  this  concept, 
binding  the  needed  hardware  to  each  individual  software 
realization  as  it  is  generated.  This  results  in  the 
duplication  of  interface  hardware  each  time  a  given  external 
variable  is  referenced.  When  the  design  program  encounters 
an  input  or  output  primitive,  it  first  generates  the 
necessary  software  followed  by  the  corresponding  support 
hardware.  The  variable  name  is  not  compared  to  previously 
generated  I/O  interfaces.  Therefore,  multiple  sensing  of  an 
input  or  issuing  of  an  output  will  cause  the  generation  of 
an  input /output  port  for  each  reference. 

2 .  Symbol  Table  Listing  of  I/O  Variables 

To  prevent  redundant  generation  of  interface  hardware 

a  record  of  each  I/O  variable  must  be  kept  which  stores  each 
variable  as  it  is  defined  by  the  user  during  the  intital 
design  specification.  This  can  be  done  very  effectively  by 
establishing  a  symbol  table  which  will  maintain  the  status  of 
each  input/output  variable  and  the  hardware  generated  for  it. 


Rather  than  producing  the  hardware  required  from  the  software 
primitives  representing  each  I/O  operation,  it  will  be 
generated  after  the  user  defines  the  external  variables  in 
the  environment  section  of  the  design  specification.  This 
feature  will  be  incorporated  in  the  input  module  being 
developed  for  the  system. 

C.  ANALOG  TO  DIGITAL  INTERFACING 

Processor  Controlled  Conversions 

The  realization  library  available  in  the  current 
implementation  of  the  design  system  allows  communications 
between  dedicated  microprocessor  systems  and  external  analog 
signals  using  a  Burr-Brown  ADC 8 2 AG  eight  bit  analog  to 
digital  converter.  Synchronization  among  contingency /task 
pairs  is  a  critical  factor  in  accurately  generating  digital 
filter  output.  Therefore,  the  sampled  data  input  must  be 
available  from  the  a  to  d  converter  on  a  periodic  basis  as 
requested  by  the  main  processor.  The  converter  available  in 
the  Intel  8080  realization  volume  is  controlled  by  an 
independent  clock  and  is  therefore  not  suitable  for  such  an 
application.  The  addition  of  an  analog  to  digital  converter 
suitable  for  digital  filtering  applications  is  therefore 
warranted . 

The  device  chosen  for  inclusion  in  the  realization 
library  is  the  Analog  Devices  AD570  and  is  an  appropriate 
selection  for  several  reasons.  The  most  significant  of  these 
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is  its  ability  to  perform  a  conversion  when  directed  by  an 
outside  signal.  Connections  are  provided  for  both  a  data 
convert  signal  and  an  associated  output  to  indicate  the 
completion  of  the  conversion  operation.  The  circuits 
required  for  interfacing  the  device  to  the  Intel  808Q  data 
bus  are  provided  in  the  manufacturer's  specification  sheets 
[Ref.  19],  The  device  is  also  capable  of  accepting  unipolar 
or  bipolar  data  input. 

2 .  Library  Entry  Development 

The  generation  of  the  library  entry  for  this  device 
is  a  simple  task.  Because  the  use  of  hardware  in  any 
realization  is  controlled  through  software  primitives,  a  new 
software  entry  was  generated  as  well.  This  new  primitive  is 
called  S.ANAIN.  User  specified  parameters  for  the  primitive 
are  determined  by  those  required  for  inclusion  of  the 
necessary  hardware  realizations.  These  consist  of  the  name 
of  the  output  signal,  the  maximum  and  minimum  values  that  it 
can  assume,  and  the  number  of  bits  in  the  digital  represen¬ 
tation.  These  parameters  represent  the  data  required  for 
correct  generation  of  the  a  to  d  converter  hardware.  The 
3dB  bandwidth  of  the  input  signal  is  also  required  in  order 
to  produce  a  buffer  amplifier  to  the  input  of  the  device. 

The  high  voltage  specification  is  used  to  determine  the  gain 
provided  by  the  buffer  amplifier.  The  interface  hardware 
listed  in  the  specification  sheets  is  not  required  for 


implementation  by  the  design  system.  The  inclusion  of  the 
primitive  S.SENSECOND  provides  a  suitable  interface  for  the 
device . 

The  a  to  d  converter  is  listed  in  the  hardware  primi¬ 
tive  H.ADC2  and  contains  the  pin  connections  necessary  to 
implement  the  device.  The  low  voltage  limit  is  evaluated  to 
determine  if  the  interface  for  bipolar  input  is  needed.  The 
title  line  for  both  of  these  primitives  is  copied  to  the 
index  of  available  realizations  that  is  at  the  beginning  of 
the  volume.  The  attributes  of  time,  storage  requirements, 
and  inclusions  are  listed  as  well.  A  complete  listing  of 
each  primitive  is  contained  in  Figure  19  and  20,  respectively. 


59 


s. ana  in  (siq*  h , 1 #b:0»  8,50, -25#  1 , 1 00 : , , , , ,  38 16, 384  7) 

com  Drimitive  to  define  orocessor  controlled  analog  input 
com  list  s  inout»hi  volt  limit»1o  volt  limit»3db  rolloff: 
com  bits#voltage  limits, bw  limits: 

com  time»stor,ext»calc,inclraddrs 

com  added  bv  m.r.  heilstedt,  may  1^83 
name  SnaOO 1 -SnaO 30 
begin  stext 

♦analog  input  channel  for  signal  <sig>» 

Grange  <h>  to  <1>  volts 
♦  3db  rolloff  at  <b>  khz 
endtext 

incl  h.conn-al  ( <Sna006> , <$na007>, gnd, <s i q> : : ) 
com  select  gain  for  buffer  amo  to  match  inout  range 


i  f 

<!> 

.It. 

0 

skip  4 

if 

<h> 

.  1  e . 

10 

skio 

6 

i  f 

<h> 

25 

skio 

8 

i  f 

<h> 

•  1  ©  • 

50 

skio 

10 

com  set  buffer  amp  if  bipolar  output 
if  <1 >  .ge.  -5  ski o  3 
if  <1>  .ge.  -12  skip  5 
if  <1>  .eq.  -25  skio  7 

com  gain  1.0  (written  101  for  inout  range  0,+10;+-5  volts 
incl  h.bufframp  ( < $na006> * <$na007> ♦ <$na005> » 1 0 , <b> : : ) 
skip  5 

com  gain  2.5  (written  25)  for  inout  range  0,+25»+-12.5  volts 
incl  h.bufframo  (<Jna006>» <$na007>, <Sna005>,25, <b>: : ) 
skip  2 

com  qain  5.0  (written  50)  for  inout  range  0,+50;+-25  volts 
incl  h.bufframp  ( < $na006> , <$na007> , < Sna005> > 50 , <b> : : ) 
incl  h.adc2  ( <Sna0 05> / <s i g> ,  <h> , < 1 > : 8 : ) 

call  s.sensecond  ( <s i g> : 8 , 1 28 ) 
com 


Figure  19.  Analog  Input  Primitive 


h ,adc2  ( i n#  out »  h  #  1 :0#8:###*#3898#3899) 

com  Dpimitive  to  define  8  bit  processor  controlled  adc 

com  list  =  i nput #output * h i  volt  limit*  lo  volt  limit: 

com  bits:lat#pwr*chips#calc»incl*addrs 

com  added  by  m.  r.  heilstedt#  may  1983 

name  $na031 -$na060 

i nc 1  h. trimoot  (gnd# <$na09 1 > * < i n> *  200 : ) 
begin  htext 

a/d  converter*  8  bit 

device  is  analog  devices  ad570#  ic  <icn> 
connec  t i ons : 
pin  1  =  n.c . 

pin  2  =  <out>C0)  (lsb) 

pin  3  -  <out>(l)  (21sb) 

pin  9  =  <out>(2)  (31sb) 

oin  5  =  <out>(3)  (91sb) 

pin  6  =  <out>(9)  (51sb) 

oin  7  =  <out>(5)  (61sb) 

oin  8  s  <out>(6)  (71sb) 

pin  9  =  <out>(7)  (msb) 

pin  10  =  +  5  volts 

oin  11  s  conv  (blank  and  .not.  convert) 

oin  12  =  -15  volts 

oin  13  =  <$na091>  (analog  input) 

pin  19  s  gnd  (analog  common) 

endtext 

if  <1 >  .It.  0  skip  6 

com  if  unipolar  input#  make  direct  connections 
begin  htext 

oin  15  =  gnd  (bipolar  offset) 

oin  16  s  gnd  (digital  common) 

endtext 
skip  11 
begin  htext 

oin  15  =  <Jna099>  (bipolar  offset) 
oin  16  s  <$na098>  (digital  common) 

endtext 

incl  h.nand9  ( 0 > 1 # *5v * <$na098> * <5na095 :  : ) 

i  nc 1  h. invert  (<$na095>, <$na096>: : ) 

incl  h.diode-sw  ( <$na096> * <$na097> : : ) 

incl  h.diode-sw  ( <$na097> * <$na098> : : ) 

incl  h.diode-sw  ( <5na098> * <$na099>  :  : ) 

incl  h.resmfgtrwt (-15v#  <5na098>#  30000: ) 

begin  htext 

oin  17  s  <jr  (data  ready) 

oin  18  =  n.c. 

endtext 

calc  icn=icntl 


Figure  20.  Analog  to  Digital  Converter  Primitive 


V.  DESIGN  EXAMPLES 


The  ability  of  the  design  system  to  produce  an  acceptable 
realization  is  demonstrated  using  a  realistic  problem  taken 
from  the  article  by  Nagle  and  Nelson  [Ref.  20].  The 
implementation  considered  is  a  fourth  order  digital  rate 
filter  with  the  transfer  function 


1.0+0.390244z“1-1.24247z“2+0 . 344333z"3+1 . 77 044 z"4 


1.0-3.02828Z  +3 . 53682z  -1.88867z  +0.397506z 


This  expression  is  redefined  as  a  combination  of  second 
order  modules  from  which  the  difference  equations  necessary 
to  generate  the  filtering  function  in  real  time  are  derived 
The  choice  of  implementation  algorithm  for  these  modules  is 
a  difficult  one  because  none  of  the  four  methods  presented 
in  chapter  three  is  clearly  superior.  Each  of  the  direct 
forms  consists  of  the  same  number  of  multiplication  and 
addition/ subtraction  operations.  However,  the  direct  form 
one  algorithm  possesses  the  minimum  storage  requirements  of 
the  four  candidates  and  will  therefore  be  used  in  the 
development  of  the  example  designs. 
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A.  CASCADE  REALIZATION 


The  cascade  implementation  of  the  proposed  filter 
consists  of  n/2  or  two  second  order  sections.  The  equivalent 
expression  is 

H(z)  =  H1(z)  .  H2(z) 

l+1.7906z~1+1.2872z~2  #  l-2.188z~1+1.3832z~2 
l+1.823z“1+0 .8959z-2  1+1. 2054z"1+0 . 4438z“2 

1 .  CSDL  Description 

Despite  the  fact  that  the  high  level  design  language 
CSDL  will  not  be  implemented  in  the  input  specification 
module,  it  remains  a  useful  tool  in  the  current  limited 
implementation  of  the  system,  providing  a  means  of  trans¬ 
lating  the  problem  statement  into  its  intermediate  form. 
Therefore,  the  first  step  in  the  development  of  the  desired 
filter  implementation  is  to  express  the  problem  in  the 
syntax  of  CSDL. 

The  identification  section  of  the  listing  contains 
all  of  the  administrative  data  associated  with  the  problem, 
such  as  the  designer’s  name,  date  of  creation,  and  project 
title.  It  is  as  follows: 

IDENTIFICATION; 

Designer:  M.  R.  Heilstedt 

Date:  4-7-83 

Project:  Sample  Cascade  Realization 


The  environment  section  is  more  easily  determined 
after  the  completion  of  the  other  subsections  in  the  CSDL 
listing.  It  will  therefore  be  developed  last. 

The  generation  of  the  contingency /task  pairs  is 
accomplished  with  the  aid  of  the  flowchart  shown  in  Figure 
21.  Each  of  the  decision  diamonds  represents  a  contingency 
to  be  tested,  while  the  process  boxes  associated  with  them 
are  their  corresponding  tasks.  The  contingency/task  pairs 
and  their  order  of  specification  in  the  contingency  list  are 
therefore  determined  by  the  sequence  of  execution  indicated 
in  the  flowchart . 

The  first  contingency /task  apir  determines  when  a 
sample  of  the  input  signal  must  be  taken.  Because  this  must 
be  accomplished  in  accordance  with  the  sampling  theorem,  the 
"dummy"  contingency  EVERY  is  used  to  force  the  generation  of 
a  sample  once  each  period  of  the  sampling  frequency. 
Arbitrarily  choosing  a  sampoing  interval  of  sixty  milliseconds 
the  first  pair  is 

EVERY  60MS  do  SAMPLE. 

The  second  contingency/task  pair  tests  for  the 
availability  of  data  from  the  analog  to  digital  converter. 
Because  the  Analog  Devices  AD570  takes  approximately  twenty- 
five  microseconds  to  perform  a  conversion,  a  test  for  the 
presence  of  data  at  its  outputs  is  required.  This  is 
accomplished  by  sensing  the  value  of  the  single  bit  input. 


Figure  21.  Flowchart  of  Cascade  Implementation 


data  ready.  The  variable  DR  is  chosen  to  represent  this 
control  line.  The  value  of  this  input  is  normally  zero,  but 
will  change  to  logical  one  while  the  conversion  process 
takes  place.  A  return  to  zero  indicates  that  the  operation 
has  been  completed.  A  sample  of  the  input  signal  is  required 
at  intervals  determined  by  the  period  of  the  sampling 
frequency.  The  second  contingency/ task  pair  is  therefore 

When  READY : 60MS  do  READ. 

The  third  contingency/task  pair  tests  to  insure  that 
the  reading  of  data  has  taken  place  and  when  true,  executes 
the  first  filtering  function  of  the  cascade.  The  time 
constraint  associated  with  this  pair  is  once  again  equal’  to 
the  period  of  the  sampling  frequency  because  data  for 
processing  is  produced  with  that  interval.  The  third 
contingency /task  pair  is  then 

When  NEWDATA: 60MS  do  FIL1. 

The  final  contingency /task  pair  in  the  list  tests  for 
the  generation  of  data  by  the  first  element  in  the  cascade. 
When  true,  the  second  filtering  function  is  executed. 

Because  the  contingency  can  be  satisfied  only  as  often  as 
data  is  produced  by  the  first  filtering  function,  its  time 
constraint  is  also  equal  to  the  sampling  interval.  The  final 
contingency/task  pair  is 

When  ONEOUT : 6 QMS  do  FIL2. 
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The  procedures  section  of  the  CSDL  listing  requires 
a  specification  of  the  high  level  code  required  to  test  the 
contingencies  and  execute  the  tasks.  The  former  shall  be 
developed  first. 

The  contingency  for  EVERY  is  a  "dummy"  specification, 
used  to  force  the  execution  of  a  specified  task  at  user 
defined  intervals.  As  such,  it  has  no  high  level  listing. 

The  contingency  READY  senses  the  value  of  the  input 
line  DR.  When  DR  is  equal  to  zero,  READY  is  set  to  one. 
Because  the  analog  to  digital  converter  will  require 
twenty-five  microseconds  to  complete  the  conversion  operation, 
a  fixed  wait  corresponding  to  this  period  is  inserted  in  the 
contingency  to  guarantee  that  sufficient  time  for  the 
conversion  to  take  place  prior  to  testing  for  its  completion. 
The  contingency  is 

Contingency  READY 

Binary  READY,1; 

Wait  2  5US ; 

Sense(DR)  ; 

If  DR=0  then  READY: =1; 

Exit  READY 

The  third  contingency,  NEWDATA,  tests  for  the  value 
of  a  flag,  called  REDDAT,  which  is  set  at  the  conclusion  of 
the  task  READ.  A  value  of  one  indicates  that  data  has  been 
read  and  is  awaiting  processing.  The  contingency  is 

Contingency  NEWDATA 

Binary  NEWDATA, 1;  f 

If  REDDAT =1  then  NEWDATA: =1; 

REDDAT :=0; 

Exit  NEWDATA 


The  final  contingency,  ONEOUT,  tests  for  the 
availability  of  data  from  the  first  filtering  function.  When 
the  value  of  the  flag  D0NE1  iw  set  to  one  at  the  end  of  task 
FIL1,  data  is  ready  for  processing  by  task  FIL2.  The  final 
contingency  is 

Contingency  ONEOUT 

Binary  oneout,!; 

If  D0NE1=1  then  ONEOUT :=1; 

D0NE1 : =0 ; 

Exit  ONEOUT 

In  addition  to  generating  the  filter  output,  the  tasks 
also  perform  internal  operations,  such  as  the  setting  of 
flags,  which  are  invisible  to  the  user.  The  first  task  to 
be  implemented  is  SAMPLE ,  which  is  used  to  produce  the  signal 
required  to  initiate  a  conversion  by  the  analog  to  digital 
converter.  This  is  done  by  setting  the  control  signal,  CONV, 
to  zero,  issuing  it  as  an  output,  resetting  its  value  to 
logical  one,  and  issuing  it  again.  The  listing  for  the  task 
is 

Task  SAMPLE 

CONV : =0  ; 

Issue( CONV )  ; 

CONV : =1 ; 

Issue (CONV)  ; 

Exit  SAMPLE 

Task  READ  senses  the  value  of  the  input  variable  as 
determined  by  the  analog  to  digital  converter.  It  must  also 
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set  the  contingency  READY  to  zero  and  the  data  available 
flag,  REDDAT,  to  one.  The  high  level  description  of  READ  is 


Task  READ 

READY :=0; 
Sense(X) ; 
REDDAT :=1; 
Exit  READ 


The  tasks  that  generate  each  of  the  filtering 
operations  are  similar  in  design  and  will  therefore  be 
developed  together.  FIL1  resets  contingency  NEWDATA  to  the 
false  condition.  The  filtering  algorithm  is  executed,  with 
the  result  being  stored  in  a  global  variable.  The  flag 
D0NE1  is  set  to  one  to  indicate  the  completion  of  the 
computation,  and  the  intermediate  variables  used  are  updated 
for  use  in  the  next  computation.  FIL2  first  sets  contingency 
ONEOUT  to  zero.  It  performs  the  same  general  calculations 
using  different  coefficients,  and  issues  the  result  of  its 
computations  as  the  system  output  variable  Y.  The  task 
listings  for  FIL1  and  FIL2  are 


Task  FIL  1 

NEWDATA: =0; 

Mil : =X-( 1 . 8  2  3*M1 2)-(0.8959*M13) ; 
Y1:=M11+(1.7906*M12)+(1. 2872*M13) ; 
D0NE1 : =1 ; 

M13 : =M12 ; 

M12 : =M11 ; 

Exit  FIL1 


Task  FIL2 

ONEOUT : =0 ; 

M21:=Y1-(1.2054*M22)-(0.4438*M23) ; 
Y:=M2 l-(2. 188 *M22)+(1. 3832 *M23); 
Issue( Y) ; 

M2  3 : =M22 ; 

M22 : =M21 ; 

Exit  FIL2 


With  the  completion  of  the  contingency  and  task 
specifications,  the  environment  section  can  be  defined.  Eac 
of  the  input  and  output  variables  must  be  listed  first.  The 
are  represented  by  X  and  DR,  and  Y  and  CONV,  respectively. 
The  arithmetic  variables  are  the  global  parameters  used 
within  the  CSDL  listing  by  multiple  contingencies  and  tasks. 
Listed  with  each  variable  is  a  description  of  the  type  of 
signal  it  is  and  the  number  of  bits  required  to  represent 
it.  The  completed  environment  section  is 


ENCIRONMENT; 

Input:  X,8,TTL;DR,1,TTL; 

Output:  Y , 8 ,TTL ; CONV ,1 ,TTL ; 

Arithmetic:  Y1 , 24 : Y2 , 24 ;M11 , 24 ;M13 , 24 ; 

M21,24:M22,24;M31,24;M32,24; 
M3 3 , 24 ; READY, 1 ;D0NE1 ,1 ; 
NEWDATA , 1 ; ONEOUT , 1 ; 

REDDAT,!; 


2 .  Intermediate  Representation 

The  generation  of  the  list  of  primitives  and  the 
IADEFL  file  for  the  example  specification  is  a  line  by  line 
translation  of  the  CSDL  listing.  The  IADEFL  file  is  derived 
from  the  identification  section  and  the  contingency  list.  I 
is  produced  by  extracting  the  names  of  the  contingency/task 


pairs  and  their  timing  constraints  from  the  former  and 
listing  them  individually.  In  addition  to  the  period  of  the 
contingency,  the  user  may  also  specify  the  maximum  time 
allowed  for  the  execution  of  the  contingency,  the  maximum 
allowed  duration  of  the  task,  the  global  order  of  the 
contingency/task  pair,  its  priority,  and  the  maximum  allowed 
time  duration  of  any  timed  block  within  the  contingency  and 
its  maximum  allowed  duration.  The  metric  upon  which  the 
design  is  to  be  generated  is  not  part  of  the  CSDL  specifi¬ 
cation  but  is  included  in  the  IADEFL  file.  The  available 
criteria  for  design  selection  are:  first  successful 
realization  produced,  realization  with  least  power  require¬ 
ments,  and  least  ccstly  realization.  Because  it  is  the  only 
one  implemented,  the  first  metric  was  chosen.  The  design 
identification  data  follows  the  listing  of  the  design 
criteria. 

The  procedures  section  of  the  CSDL  listing  contains 
the  information  needed  to  generate  the  primitive  list 
specification  of  the  design  problem.  Each  line  of  the  high 
level  expression  is  translated  into  one  or  more  lines  of 
intermediate  code.  In  order  to  manually  generate  the 
primitive  problem  representation,  the  user  must  have 
available  a  listing  of  the  index  to  the  realization  volume 
used.  Each  of  the  entries  in  the  intermediate  specification 
is  derived  from  the  title  lines  of  the  available  primitives. 
The  process  of  translating  the  CSDL  listing  into  its 
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intermediate  form  is  tedious  and  repetitious  and  will  not  be 
detailed  here.  However,  a  short  example  is  instructive.  The 
line  to  be  translated  is  a  mathematical  expression  taken 
from  the  first  filtering  task,  FIL1.  The  example  statement 
is 

Mil : =X-( 1 . 8  23*M12 ) -( 0 . 8959*M13 ) 


The  corresponding  primitive  code  is  produced  by  working  from 
right  to  left  in  this  statement.  Each  of  the  expressions  in 
parentheses  is  derived  first.  The  results  of  these  operations 
are  then  combined  as  indicated  by  the  remaining  operations  to 
attain  the  desired  result.  The  constants  must  also  be 
specified  as  variables  within  the  listing.  Each  variable 
name  used  requires  an  additional  statement  for  the  genera¬ 
tion  of  a  storage  location.  The  primitive  list  for  this 
equation  is 


s . fmul 

(mla ,all ,ml2 : 24 , 24 , 24 ) 

s . fmul 

(mlb , al 2, ml 3:24, 24, 24) 

s . fsub 

(ma  ,mla  ,mlb :  24 , 2 4 , 24- ) 

s . float 

( sigin ,insig : 8  ) 

s .fsub 

(mil ,sigin,ma:24,24,24) 

s .  var 

(mla : 24 ) 

s  .var 

(ml2: 24) 

s  .var 

(mlb: 24) 

s  .var 

( ml 3:24) 

s .  var 

(ma : 24 ) 

s  .var 

( insig : 24 ) 

s . fcons 

(all, 0,0, 0,0:24) 

s .fcons 

(al2 ,0,0,0,0:24) 

The  IADEFL  file,  primitive  list,  CSDL  problem  statement,  and 
design  system  output  are  includes  in  Appendix  A. 


B.  PARALLEL  REALIZATION  ' 

The  parallel  realization  of  the  fourth  order  digital 
rate  filter  also  consists  of  two  second  order  modules.  The 
filtering  function  to  be  generated  is  found  by  taking  the 
partial  fraction  expansion  of  the  original  transfer  function 
The  resulting  expression  is 

H(z)  =  1.0  -  S.O1--26402  1-^-7475z  1 

1-1.8230Z  +0.8959z 

4  q1.6734z~1-1. 5692z~2 
l-1.2053z_1+0.4437z“2 

It  is  illustrated  by  the  flowchart  of  Figure  22. 

1.  CSDL  Description 

The  CSDL  description  of  the  parallel  realization  is 
similar  to  its  cascade  counterpart.  Therefore,  only  those 
areas  in  which  they  differ  are  presented  in  detail.  The 
identification  section  for  the  parallel  realization  is 

IDENTIFICATION; 

Designer:  M.  R.  Heilstedt 

Project:  Sample  Parallel  Filter  Problem 

Date:  7  April  1983 

The  dummy  contingency  EVERY  is  once  again  used  to 


force  the  sampling  of  the  input  signal  according  to  the 
interval  determined  by  the  sampling  frequency.  The  second 
contingency  tests  for  the  completion  of  the  conversion  and 


START 


Figure  22.  Flowchart  of  Parallel  Implementation 


executes  the  task  which  reads  the  digital  data  input.  Each 
of  these  contingency/task  pairs  is  identical  to  those  of  the 
same  name  in  the  cascade  problem  specification. 

The  remaining  entries  in  the  contingency  list  are 
unique  to  the  parallel  problem  and  are  used  to  implement 
the  FORK  construct.  The  first  of  these  tests  for  the 
availability  of  data  for  processing.  A  true  condition 
executes  the  task  FORK,  which  sets  two  variable  flags,  VAR1 
and  VAR2 ,  to  one.  Contingencies  Cl  and  C2  test  the  values 
of  these  flags.  When  true,  Cl  permits  the  execution  of  FIL1 
the  first  of  the  parallel  filtering  functions.  Similarly, 
the  second  filtering  slgorithm,  FIL2,  is  executed  when  the 
contingency  C2  is  satisfied.  FIL1  and  FIL2  each  place,  their 
data  in  individual  queues  and  increment  the  value  of  their 
respective  queue  counters.  Contingency  C2  tests  for  the 
existence  of  each  queue  and  if  both  are  present,  causes  the 
execution  of  task  PLUS.  Once  again  choosing  an  arbitrary 
sampling  interval  of  sixty  milliseconds,  the  contingency 
list  is 

CONTINGENCY  LIST; 

EVERY  60  MS  do  SAMPLE 
When  DATA : 6  QMS  do  READ 
When  LOOP : 60MS  do  FORK 
When  Cl  :60MS  do  FIL1 
When  C2  :60MS  do  FIL2 
When  C3  :60MS  do  PLUS 

Using  the  terminology  introduced  in  chapter  four,  FIL1  and 
FIL2  represent  the  forked  tasks  and  PLUS  is  the  joined  task. 


The  algorithm  for  contingency  DATA  is  identical  to 
that  used  for  its  counterpart  in  the  cascade  implementation. 
Contingency  LOOP  tests  the  value  of  variable  GO  which  is  set 
to  one  by  the  execution  of  task  READ.  When  the  contingency 
is  satisfied,  it  permits  task  FORK  to  be  executed.  The 
contingency  is 


Contingency  LOOP 

Binary  L00P,1; 

If  G0=1  then  L00P:=1; 
G0  =  0; 

Exit  LOOP 


Contingencies  Cl  and  C2  perform  similar  functions, 
using  the  variables  VAR1  and  VAR2,  respectively.  When  true, 
Cl  permits  execution  of  FIL1  and  C2  permits  execution  of 
FIL2.  Both  of  the  test  variables,  VAR1  and  VAR 2 ,  are  set  by 
execution  of  task  FORK.  The  contingencies  are 


Contingency  Cl 

Binary  Cl,l; 

If  VAR1=1  then  Cl=l; 
VAR1  =  0 ; 

Exit  Cl 

Contingency  C2 

Binary  C2,l; 

If  VAR2=1  then  C2=l; 
VAR2  =  0  ; 

Exit  C2 


Contingency  C3  tests  the  value  of  the  queue  counters 
to  determine  if  data  generated  by  tasks  FIL1  and  FIL2  is 
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waiting  to  be  processed.  If  both  queues  contain  data 
entries,  task  PLUS  is  executed.  The  CSDL  listing  for  C3  is 

Contingency  C3 

Binary  C3,l; 

If  CTRl.gt.O  .and.  CTR2.gt.O  then  C3:=l; 

Exit  C3 

The  tasks  SAMPLE  and  READ  are  identical  to  the  tasks 
already  described  for  the  parallel  filter  implementation 
problem. 

Task  FORK  sets  the  flags  tested  by  contingencies  Cl 
and  C2  to  one  to  represent  the  true  condition.  It  also 
resets  the  value  of  its  associated  contingency,  LOOP,  to 
zero.  The  task  listing  is 

Task  FORK 

LOOP:=0; 

VAR1 : =1 ; 

VAR 2 : =1 ; 

Exit  FORK 

It  is  important  to  note  that  the  contingency/task  pairs  that 
follow  cannot  be  satisfied  if  the  LOOP/FORK  pair  is  not  true 
as  well.  Because  the  tasks  that  perform  the  filtering 
operations  are  included  in  this  group,  FORK  effectively 
controls  the  execution  of  the  filtering  function. 

The  CSDL  descriptions  of  tasks  FIL1  and  FIL2  are 
similar  to  those  used  in  the  cascade  realization,  "’he  direct 
form  one  algorithm  is  again  used  to  implement  the  filtering 
function.  The  results  produced  by  each  of  the  tasks  are 


placed  in  individual  data  queues  and  their  associated 
counters  are  incremented.  The  high  level  descriptions  of 
the  tasks  are 


Task  FIL1 

Cl : =0 ; 

Mil : =X+ (  1. 8  23*M12  ) - (  0 . 8  96*M13 ) ; 
Q1B :  =  ( 1 . 264*M11 ) - ( 1 . 74 8*M12 )  ; 
M13 : =M12 ; 

Ml 2 : =M11 ; 

CTR1 : =CTR1+1 ; 

Exit  FIL1 

Task  FIL2 

C2 :  =  0  ; 

M21:=X+(1.205*M22)-(0.444*M23) ; 
Q2B:=(1.673*M21)-(1.569*M22) ; 
M23:=M22; 

M22 : =M21 ; 

CTR2 : =CTR2+1 ; 

Exit  FIL2 


Task  PLUS  adds  the  data  generated  by  FIL1  and  FIL2  to 
produce  the  output  of  the  device.  The  queue  counters  are 
decremented  following  the  output  of  the  filtered  signal.  The 
CSDL  description  of  PLUS  is 


Task  PLUS 

C  3  :  =  0  ; 

Q1A : =Q1B ; 

Q2A: =Q2B ; 

Y: =1 . 0-( 8 . 0*Q1A) -( 4 . 0*Q2A) ; 
Issue(Y) ; 

CTR1 : =CTR1-1 ; 

CTR2 : =CTR2-1 ; 

Exit  PLUS 


Now  that  the  contingencies  and  tasks  for  the  parallel 


problem  have  been  defined,  the  environment  section  entries 
can  be  determined. 


ENVIRONMENT; 

Input:  X,8 ,TTL;DR,1 ,TTL; 

Output:  Y , 8 ,TTL ;C0NV ,1 ,TTL ; 

Arithmetic:  Q1A, 24 ;Q1B , 24 ;Q2A, 24 ;Q2B , 24 ; 

Mil , 24 ;M12 , 24 ;M13 , 24 ;M21 , 24 ; 
M22,24;M23,24;G0,1  ;VAR1 , 8  ; 
VAR 2, 8; LOOP, 8; DATA , 8 ; 


2 .  Intermediate  Representation 

The  IADEFL  file  and  the  primitive  list  are  generated 
from  the  CSDL  listing  as  detailed  in  the  development  of  the 
input  specification  for  the  cascade  problem.  Appendix  B 
contains  these  listings,  as  well  as  the  complete  CSDL 
specification  of  the  problem  and  the  output  date  generated 
by  the  design  program. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  feasibility  of  automating  the  production  of 
microprocessor-based  digital  filters  has  been  demonstrated. 

The  use  of  applications  problems  has  been  shown  to  be  a 
logical  and  productive  step  in  the  development  of  the  design 
system.  The  modifications  and  additions  made  as  a  result  of 
this  research  have  increased  the  versatility  of  the  current 
version  of  the  design  program  and  identified  areas  for 
future  study. 

The  value  of  the  design  system  described  and  tested  in 
this  thesis  is  readily  apparent.  The  time  now  required  to 
manually  produce  the  hardware  and  software  necessary  to 
support  dedicated  microprocessor  systems  can  be  more 
productively  spent  if  the  computer  can  be  programmed  to 
accomplish  this  for  us.  Further  development  of  the  design 
system  is  therefore  warranted. 

B.  RECOMMENDATIONS 

1 .  Implementation  Language 

The  current  version  of  the  design  program  is  written 
in  FORTRAN.  This  language  places  rigid  requirements  on  the 
format  of  the  input  data  files  used  by  the  system  and 
restricts  the  organization  of  the  program  code.  Modularization 


of  the  program  algorithm  into  specific  levels  of  operation  is 
not  possible.  The  use  of  a  language  such  as  Pascal,  PL/1, 
or  Ada  would  permit  this  type  of  program  organization, 
allowing  additions  and  modifications  to  be  made  more  easily. 
Validation,  verification,  and  testing  of  the  design  program 
would  also  be  enhanced. 

2 .  Validation  of  Current  Program 

The  design  program  installed  on  the  Digital  Equipment 
VAX  11/780  was  received  on  magnetic  tape  from  Lawrence 
Livermore  Laboratory.  Because  of  a  compatibility  problem 
between  the  machine  that  produced  the  tape  and  that  which 
read  it,  sporadic  errors  occurred  in  the  copying  of  the  tape 
onto  the  VAX.  These  errors  consisted  primarily  of  incorrect 
branch  statements ,  but  erroneous  variable  references  were 
also  found.  These  problems  were  of  sufficient  severity  to 
prevent  realizations  from  being  generated  when  in  fact  they 
were  possible.  Therefore,  a  considerable  amount  of  time 
was  spent  debugging  the  design  program  and  in  fact  became  a 
major  portion  of  the  effort  to  produce  this  thesis.  A 
thorough  validation  of  the  program  remains  to  be  accomplished 
That  which  has  been  done  thus  far  has  been  performed  to 
identify  the  source  of  observed  errors  in  program  execution. 
Due  to  the  nature  of  the  inconsistencies  found  it  is 


reasonable  to  assume  that  those  subroutines  not  checked 
contain  errors  as  well.  A  line  by  line  manual  comparison  of 


the  listing  of  the  design  program  installed  on  the  VAX  11/780 
with  a  copy  of  the  code  that  is  known  to  be  correct  is  the 
best  method  for  determining  the  correctness  of  the  current 
imp  1  erne  nt  at  ion . 

3 .  Realization  Library 

The  choice  of  implementation  language  for  the  design 
program  will  also  determine  the  format  of  the  entries  in  the 
realization  volume.  Its  current  organization  makes  additions 
and  modifications  a  difficult  task.  Independent  of  the 
implementation  language  and  library  format  is  the  need  for 
an  algorithm  to  allow  library  updates  by  the  user.  This  type 
of  program  would  save  considerable  time  in  the  development  of 
the  hardware /software  database. 

The  organization  of  the  library  merits  further  study. 

In  keeping  with  the  concept  of  hardware  binding,  the  user 
specifies  his  problem  in  terms  of  software  primitives.  The 
support  hardware  necessary  is  automatically  generated. 

Because  the  design  algorithm  is  intended  to  search  only  one 
volume  at  a  time  for  the  needed  primitives,  duplication  of 
hardware  between  volumes  can  occur.  To  prevent  this  redundancy 
the  creation  of  a  global  realization  volume  of  commonly  used 
hardware  primitives  would  be  helpful.  After  an  unsuccessful 
search  of  the  current  microprocessor  volume,  the  program 
would  scan  the  global  listing  in  an  attempt  to  locate  the 
needed  primitive.  Failure  to  produce  a  realization  would 


occur  only  after  an  unsuccessful  search  of  both  the  current 


and  global  volumes.  The  content  of  the  global  volume  would 
not  be  limited  to  modular  hardware.  Individual  components, 
wuch  as  resistors  and  capacitors,  could  be  included  as  well. 
The  principle  of  hardware  binding  is  not  violated  as  long  as 
user  access  to  these  entries  is  only  allowed  through  software 
primitives . 

The  ultimate  goal  of  the  design  system  is  to  produce 
a  physical  hardware  realization  of  the  problem  specification. 
In  such  a  system,  the  current  program  would  be  the  first 
module  in  a  group  of  three  or  more  that  would  generate  the 
layout  for  the  device  program  read  only  memory,  with  the 
monitor  and  design  program,  and  assemble  the  hardware.  The 
program  code  required  to  implement  such  a  system  would 
necessitate  the  listing  of  the  software  primitives  in  terms 
of  microprocessor  operational  codes  rather  than  assembly 
language.  Therefore,  the  organization  and  incorporation  of 
op  code  based  volumes  in  the  realization  library  must  be 
investigated . 

4 .  Interrupt  Driven  Monitors 

The  theory  for  the  implementation  of  an  interrupt- 
driven  monitor  has  been  developed  but  inclusion  in  the 
design  system  remains  to  be  accomplished.  The  availability 


of  such  a  strategy  as  a  user  selectable  alternative  will 
greatly  increase  the  potential  of  the  design  system. 


More  applications  problems  are  required  to  determine 
the  remaining  shortcomings  in  the  design  system.  Signal 
processing  implementations  involving  modulation/demodulation 
techniques,  as  well  as  high  speed  special  purpose  hardware 
realizations,  remain  to  be  tested.  The  availability  of  the 
interrupt-driven  monitor  strategy  can  be  expected  to  create 
additional  test  applications. 

6 .  Documentation 

The  documentation  currently  available  is  inadequate 
to  permit  rapid  familiarity  with  the  program.  Therefore, 
reference  data  must  be  maintained  which  details  the 
development  and  implementation  of  the  system.  This  will  aid 
subsequent  research  and  provide  the  basis  for  a  user's 
manual  when  a  commercially  use 'V  e  design  system  is  produced 


APPENDIX  A 


DATA,  CASCADE  REALIZATION 


This  appendix  contains  the  input  and  output  files  for 
the  cascade  realization  developed  in  chapter  five. 


A.  INPUT  DATA 

1.  CSDL  Listing 

"This  is  the  CSDL  description  of  a" 

"fourth  order  digital  'rate'  filter.  It  is" 

"taken  from  the  article  by  Nagle  and" 

"Nelson  which  appeared  in  the  February  1981  issue" 
"of  COMPUTER  magazine.  It  is  implemented  as" 

"the  cascade  of  two  second  order  sections,  which" 
"are  expressed  using  the  direct  form  one  algorithm." 

IDENTIFICATION; 

Designer:M.  R.  Heilstedt 
Date : 4-7-83 

Project:  Sample  Cascade  Filter  Problem 

ENVIRONMENT; 

Input : X , 8 ,TTL;DR , 1 ,TTL ; 

Output :Y, 8, TTL;C0NV,1,TTL; 

Arithmetic : Y1 , 8 ; Y2 , 8 ;M11 , 8 ;M13 , 8 ;M21 , 8 ;M22 , 8 ; 
M31,8;M32,8;M33,8 ; READY ,1 ; 

D0NE1 ,1 ; 

CONTINGENCY  LIST; 

"Sample  signal  every  60  milliseconds" 

EVERY  60ms  do  SAMPLE 
"Check  for  conversion  completed" 

When  READY :60ms  do  READ 
"Do  first  second  order  filtering" 

When  NEWDATA :60ms  do  FIL1 
"Then  do  second  filtering" 

When  ONEOUT :60ms  do  FIL2 
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UNCLASSIFIED 


F/G  9/2 


NL 


PROCEDURES; 

Contingency  READY 
Binary  READY,1; 

Wait  25US; 

Sense(DR) ; 

If  DR=  0  then  READY: =1  fi; 

Exit  READY 
Contingency  NEWDATA 
Binary  NEWDATA, 1; 

If  REDDAT=1  then  NEWDATA :=1 
REDDAT : =0 ; 

Exit  NEWDATA 

Contingency  ONEOUT 
Binary  ONEOUT, 1; 

If  D0NE1=1  then  ONEOUT :=1; 

DONE1 :  =  0 ; 

Exit  ONEOUT 

Task  SAMPLE 
CONV : =0 ; 

Issue ( CONV )  ; 

CONV : =1; 

Issue(CONV) ; 

Exit  SAMPLE 

Task  READ 

READY: =0; 

Sense(X) ; 

REDDAT: =1; 

Exit  READ 

Task  FIL1 

NEWDATA: =0; 

Mil : =X-( 1 . 8  23*M1 2)  —  (0.89 5  9*M13 ) ; 
Y1 : =M11+ ( 1 . 7  9  0  6*M12  )+ (  1 . 28  72*M1 3 ) 
D0NE1 : =1 ; 

M13 : =M12 ; 

M12 : =M11 ; 

Exit  FIL1 

Task  FIL2 

ONEOUT :=0; 

M21:=Y1-(1.2054*M22)-(0.4438*M23) 
Y:=M21-( 2.188*M22)+(1.3832*M23) ; 
Issue( Y) ; 

M23 : =M22 ; 

M22 : =M21 ; 

Exit  FIL2 


001  isystem  : 


2 .  IADEFL  File 
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3.  Primitive  Listin; 


o 

a 

o 

o 

O 

P 

o 

o 

o 

o 

o 

o 

o 

o 

o 

p 

o 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

p 

o 

D 

p 

p 

p 

o 

o 

o 

o 

o 

D 

O 

D 

O 


********* 


********* 


generated  for:svstem 
main  ( : : ) 

start  ( : : ) 

generated  forteach 
every  (each::) 

var  (each:8) 

generated  for:data 
proc  (data::) 

f i xedwai t  (25: : ) 
sensecond  (dr:l,l) 
eg  Ot2,dr,3cOOl :8/ 1 , 1 ) 

jmpf  (9t2#31 001:8) 

ass igncons( data# 1:1*1) 


********* 


1  oc 

ex i tproc 

cons 

var 

var 

var 


) 

1*1) 


********* 


(31 001 : ) 

(data»0: 

(3c 00 1  *  0: 

(dr: 1 ) 

(data:  1 ) 

(3t2:8) 

generated  for:newdata 
proc  ( newdata : : ) 

eq  (it  3# reddat #3c002:8*  1 ,  1 ) 

jmpf  (3t3/31002:8) 

ass i gncons (reddat, 0:1,1) 
ass i gncons (newdata, 1:1,1) 

1 oc  01002:8) 

exitproc  ( newdat a ,  0  : : ) 

0c002,0:  1  ,  1  ) 

(newdata :  1 ) 

( reddat :  1 ) 

Ot3:8) 

generated  for:oneout  ********* 
proc  (oneout::) 

ea  Ot  ^ ,  done  1 , 3c  00 3 : 8 , 1  ,  1  ) 

jm of  Ota, 31 003: 8) 

ass i gncons (oneout ,1:1,1) 
assigncons( done  1 ,0:1,1) 
l  oc  01003:) 

exitproc  (oneout, 0::) 

(3c003, 1  :  1  ,  1  ) 

(oneout : 1 ) 

(donel :  1 ) 

Ot«:8) 

generated  for:sample 
proc  (sample::) 

ass i gncons (con v  ,0:1,1) 


cons 

var 

var 

var 


cons 

var 

var 

var 


********* 


s.issuecond  (conv:l,2) 
s.assigncons(conv, 1:1,1) 
s.issuecond  (conv:l,2) 
s.exitproc  (sampl e,each: : ) 

s. var  (convJl) 

t.  generated  forjpead  ****** 

s.proc  (read:) 

s.assigncons(data,0: 1,1) 
s.anain  (x, -10, 10,5:8: ) 

s.ass i qncons ( reddat ,1:1,1) 
s.exitproc  ( read , da t a : : ) 

s. var  (x:8) 

t.  generated  for:fill  ****** 

s.oroc  (fill::) 

s.assigncons(newdata,0: 1,1) 
s.fmul  (mlb,b2,ml 3:2a,2«,2«) 

s.fmul  (mla,bl,m!2:24,24,24) 

s.fsub  (ma,mla,mlb:24,24,24) 

s.  f 1  oat  ( i x, x : 8) 

s.fsub  (ml 1 , i x,ma:24,24,24) 

s.fmul  (ylb,a2,ml3:24,24,24) 

s.fmul  (vla,al , ml  2 : 24 , 2a , 24 ) 

s . f add  (ya,yla,y lb:24,24,24) 

s. f add  (vl,mll,ya:24,24,24) 

s. ass igncons( done  1 , 1:1,1) 
s.fassign  ( m 1 3, m 1 2 : 24 , 24) 

s.  f assign  (ml 2,m 1 1 : 24 ,24) 

s.exitproc  ( f i 1 1 , newdat a : : ) 
s.var  ( m 1 b : 24 ) 

s.var  ( m 1 3 : 24 ) 

s.var  (mla:24) 

s.var  (ml2:24) 

s.var  (ma:24) 

s.var  (ml  1 :24) 

s.var  ( i x :24) 

s.var  (y lb:24) 

s.var  (yla:24) 

s.var  (ya:24) 

s . f cons  (b2,0, 168, 1 14,64:24) 

s . f cons  (bl ,0, 1 70, 1 14,64:24) 

s . f cons  (a2,0, 142, 164,64:24) 

s .  f cons  (al ,0,22,228,64:24) 

t.  generated  for:fil2  ******* 

s.proc  ( f i 1 2: : ) 

s.assigncons(oneout,0: 1,1) 

S.fmul  (m 2b, m23, 622:24,24,24) 

s.fmul  (m2a,m22,b21:24,24,24) 

s.fsub  (m2a,m2a,m2b:24,24,24) 

s.fsub  (m2l,yl,m2a*24,24,24) 

s.fmul  (yb,a22,m23:24,24,24) 


s.fmul  ( yh , a2 1 , m22 : 24 , 24 , 24 ) 

s.fadd  (yh,yh,yb:24,24,24) 

s. f Sub  (y/m21,yh:2a,2«/2a) 

S  .  f  i  x  ( i  y ,  y  :  8 ) 

s.anaout  ( i y / - 1 0 / 1 0 : 8 ) 

s.fassiqn  (n23,m 22:24,24) 

s.fassign  ( *22 , m21 t 24 / 24 ) 

s.exltoroe  ( f i 1 2,oneout : : ) 

s.var  (>n2b:24) 

s.var  (m23:24) 

s.var  (m2a:24) 

s.var  (m22:24) 

s.var  ( yb : 24 ) 

s.var  ( y  h  :  24 ) 

s.var  (y : 24 } 

s.var  ( i y:8) 

s.fcons  ( b22 , 0 , 168 , 56 » 64 : 24  ) 

s.fcons  (b2 1 , 0 , 36 , 205 » 64 : 24 ) 

s.fcons  Ca22, 0,8, 49, 64:24) 

s.fcons  (a21, 0,4, 198,192:24) 

s .end  ( : : ) 


OUTPUT  DATA 

1.  Software  Listin 
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APPENDIX  B 


DATA,  PARALLEL  REALIZATION 


This  is  the  input  and  output  data  for  the  parallel 
filter  problem  developed  in  chapter  five.  The  timing 
constraints  were  increased  to  ninety  milliseconds  to  permit 
successful  generation  of  a  realization.  The  design  program 
attempted  to  generate  a  dual  processor  realization  for  the 
original  specification,  but  a  system  error  prevented 
completion  of  the  output. 


A.  INPUT  DATA 

1 .  CSDL  Listing 

"This  is  a  parallel  realization" 

"of  the  fourth  order  rate  filter" 

“as  presented  in  the  article" 

"by  Nagle  and  Nelson  in  the" 

’’February  1981  issued  of  IEEE  COMPUTER." 

IDENTIFICATION: 

Designer:  M.  R.  Heilstedt 
Date:  7  April  1983 

Project:  Sample  Parallel  Filter  Problem 

ENVIRONMENT ; 

Input : X , 8 ,TTL ;DR ,1 ,TTL ; 

Output :Y , 8 ,TTL ;CONV ,1 ,TTL ; 

Arithmetic :Y1 , 8 ;M11 , 8 ;M12 , 8 ;M13 , 8 ; 

Q 1A, 24 ;Q2A, 24 ;Q1B, 24;Q2B ,24 ;M21 ,8 ;M22 ,8 ;M23 , 8; 

GO , 1 ; VAR1 , 1 ; VAR2 , 1 ; CTR1 , 1 ; CTR2 , 1 ; 

CONTINGENCY  LIST: 

"Initiate  sequence  every  60  milliseconds: 
EVERY  60ms  do  SETFLG 
"Sample  Signal  once  each  period" 

EVERY  60ms  do  SAMPLE 


"Check  for  data  available  after  sample  converted 
When  DATA :60ms  do  READ 

"Set  flags  for  parallel  operations  when  ready" 
When  LOOP  :60ms  do  FORK 
"First  parallel  filter  function" 

When  Cl :60ms  do  FIL1 
"Second  parallel  filter  function" 

When  C2 :60ms  do  FIL2 

"Add  results  of  each  function  and  issue  sum" 

When  C3 :60ms  do  PLUS 

PROCEDURES; 

Contingency  DATA 
Binary  DATA,1; 

Sense  ( DR) ; 

If  DR=0  then  DATA : = 1 ; 

Exit  DATA 

Contingency  LOOP 
Binary  L00P,1; 

If  G0=1  then  L00P:=1; 

G0  =  0 ; 

Exit  LOOP 

Contingency  Cl 

Binaray  Cl,l; 

If  VAR1=1  then  Cl:=l; 

VAR1 : =0 ; 

Exit  Cl 

Contingency  C2 
Binary  C2,l; 

If  VAR2=1  then  C2 . =1; 

VAR2 : =0 ; 

Exit  C2 

Function  C 3 

Binary  C3,l; 

If  CTRl.gt.O  .and.  CTR2.gt.0  then  C3:=l; 

Exit  C3 

Task  SAMPLE 
CONV :  =  0  ; 

Issue(CONV) ; 

CONV : =1 ; 

Issue(CONV ) ; 

Exit  SAMPLE 


Task  READ 

DATA :  =  0 ; 

Sense (X) ; 

GO : =1 ; 

Exit  FORK 

Task  FIL1 
Cl :  =  0 • 

Mil : =X+ (1 . 82  3*M12 )-( 0 . 89  6*M13 ) 
Q1B : = ( 1. 264*Mll)-( 1 . 748*M12 ) ; 
M13 : =M12 ; 

Ml 2 : =M11 ; 

CTR1 : =CTR1+1 ; 

Exit  FIL1 

Task  FIL2 
C2L=  0 ; 

M21 : =X+ ( 1 . 205*M22 )-( 0 . 444*M23 ) 
Q2B:=(1.673*M21)-(1.569*M22) ; 
M23 : =M22 ; 

M22:=M21; 

CTR2 : =CTR2+1 ; 

Exit  FIL2 

Task  PLUS 
C3 :  =0 ; 

Q1A : =Q1B ; 

Y:=1.0-(8. 0*Q1A) - ( 4 . 0*Q2A) ; 
Issue(Y) ; 

CTR1 : =CTR1-1 ; 

CTR2 : =CTR2-1 ; 

Exit  PLUS 


Primitive  Listin 


generated  for:system  ******* 

main  (  : :  ) 

start  ( : : ) 

generated  forteach  ******* 

(each : : ) 
var  (each:8) 

generated  for:data  ******* 

oroc  (data::) 

f i xedwa it  ( 25 : : ) 

sensecond  (dr:l,l) 

eq  Ot  2 ,  dr  ,  3c  00  1 :  8 , 1  *  1 ) 

jmpf  Ot2, 31001:8) 

assigncons(data,  1 :  1 ,  1 ) 
loc  (31001::) 

exitproc  (data>0::) 

cons  (3c001,0:l) 

var  Ot2: 8) 

var  (data:l) 

var  (dr:l) 

generated  for:loop  ******* 

proc  (loop::) 

eq  (3t3,go,3c002:8, 1 , 1 ) 

jmpf  Ot3, 31 002:8) 

assi gncons ( 1 ooo# 1 : 1 t  1 ) 
assi gncons (go#  0 : 1 ,  1 ) 
loc  (31002::) 

exitproc  (loop»0::) 

cons  (9c002»  1:1) 

var  (3t3:8) 

var  (1 oop : 1 ) 

var  (go: 1 ) 

generated  for:cl  ******* 

oroc  (cl::) 

eq  (3t 4 t var 1 > 3c00 3 : 8 t  l f 1 ) 

jmpf  Ota, 31 003:8) 

assi gncons (c  1 , 1 : 1 , 1 ) 
assi gncons ( var 1 ,0:1,1) 
loc  01003::) 

ex i toroc  (c  1 , 0 : : ) 

cons  (3c 003, 1:1) 

var  Ot  a :  8) 

var  (c  1 :  l ) 

var  (var 1 : 1 ) 

generated  for:c2  ******* 

proc  (c2::) 

eq  Ot 5 ,  var 2 , 3c 00a  :  8 ,  1 , 1 ) 

jmof  (3t5,31 00a:8) 

assi gncons (c2, 1:1,1) 


********************* 


********************* 


140 


assi gncons (var2, 0:1,1) 
loc  (91004::) 

exi toroc  (c2, 0: : ) 
cons  ( 9c  004 , 1:1,1) 

var  (3t5:8) 

var  (var2:l) 

var  (c2: 1 ) 

generated  for:c3  ******* 

oroc  (c3::) 

gt  (  3t  6,  c  t  r  1 , 3c  005  :  8 , 8 , 8  ) 

jmpf  (9t 6,91 005:8) 

qt  ( 9t 6, c t r2 , 3c 005 : 8 , 8 , 8  ) 

jmpf  (9t6, 31 005:8) 

assiqnc on s(c3, 1:1,1) 


var 

cons 


sub  (c t r 1 , ct r 1 ,dec : 8 , 8 , 8) 

sub  (ctr2,ctr2,dec:8,8,8) 

loc  (31005::) 

exi toroc  (c3,0: : ) 

cons  (3c005,0:8) 

var  ( 3tb: 8) 

var  (c3: 1 ) 

var  (ctrl :8) 

var  (ctr2:8) 

cons  (dec, 1:8:) 

generated  for:samp1e  ******* 

proc  (sample::) 

ass i gncons (con v, 0:  1,1) 
issuecond  (conv:l,2) 
assi gncons (conv, 1:1,1) 
issuecond  (conv:l,2) 
exitproc  ( samp  1 e, each :: ) 
var  (conv : 1 ) 

generated  for:read  ******* 

oroc  (read::) 

assi gnc on s(data, 0:1,1) 
anain  ( i ns i g, 1 0 , - 1 0 , 1 : 8 : ) 

assi gncons ( go,  1  : 1 ,  1  ) 
exitproc  ( read , dat a  : : ) 
var  (insig:8) 

generated  for:fork  ******* 

oroc  ( f or k : : ) 

assi gncons ( 1 ooo, 0:1,1) 
ass i gncons ( var  1 , 1 : 1 , 1 ) 
assi gncons (var2, 1:1,1) 
exitproc  ( fork , 1 ooo :: ) 

generated  for:fill  ******* 

proc  (fill::) 

assi gncons ( c 1 , 0 : 1 ,  1  ) 
fmul  (ml  a, al 1 , m 1 2 : 24 , 24 , 24 ) 

fmul  ( m 1 b , a  1 2 , m 1 3 : 24 , 24 , 24 ) 


a  o  o  o  o  q  o 
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S . f SUD 

(ma#mla»mlb:24,24#24) 

s . f 1  oat 

(sigin, i ns i o: B) 

s . f add 

(ml l#sigin«ma:24r24,24) 

s. fmul 

(mla*ml2#bli :24,24,24) 

s  *  f mg  1 

( m 1 b  ?  m 1 1 , b 1 2 : 24 , 24 , 24 ) 

s. f sub 

(qlb,mlb,mla:24,24,24) 

s.fassign 

(ml 3,ml2:24,24) 

s.fassign 

(ml2,mll:24,24) 

s.add 

(ctrl#ctrl#inct0»8,8) 

s .ex i toroe 

(fill, cl::) 

s .cons 

( i nc» 1 :8) 

s.var 

( qlb : 24 ) 

s .  var 

( m 1  a : 24 ) 

s.var 

(ml2:24) 

s.var 

( m 1 b : 24 ) 

s.var 

(m  1  3: 24 ) 

s.var 

( m  a :  2  4 ) 

s.var 

(s»gin:8) 

s . f cons 

(all, 0,246, 208, 64:24) 

s . f cons 

(al2, 0,212, 239, 192:24) 

s. f cons 

(bl 1,0, 172,244, 192:24) 

s . f cons 

(bl2,0, 172, 1 14,64:24) 

t.  generated 

for: f i 1 2  ********************* 

s.oroc 

(f i 12::) 

s.ass»gncons(c2,0: 1 , 1 ) 

s. fmul 

(m2a,m22,a21:24,24,24) 

s. fmul 

(m2b,m23,a22:24,24,24) 

s . f sub 

(mb,m2a»m2b:24,24,24) 

s . f 1  oat 

(intx,insig:8) 

s . f add 

(m21,intx,mb:24,24,24) 

s . f mu  1 

(m2a,m21 ,b2l :24,24,24) 

s • f mu  1 

(m2b,m22,b22:24,24,24) 

s. f sub 

(q2b,m2a,m2b:24,24,24) 

s. f assi gn 

(m23,m22:24,24) 

s.fassign 

( m22, m2 1:24, 24) 

s .  add 

(ctr2,ctr2,inc:8,8,8) 

s.*xi tproc 

(f i 1 2,c2: : ) 

s.var 

(q2b:24) 

s.var 

(m2a : 24 ) 

s.var 

(m22 : 24 ) 

s.var 

(m2b:24) 

s.var 

(m23:24) 

s.var 

(mb : 24 ) 

s.var 

( i  nt  x : 8 ) 

s.var 

( m2 1 : 24 ) 

s • f cons 

(a21, 0,24, 235, 64:24) 

s . f cons 

(a22,0, 104,228, 192:24) 

s . f cons 

(621,0,34,205, 192:24) 

s . f cons 

(622,0,202,56,64:24) 

t.  generated 

for:p1us  ********************* 
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0 

s .proc 

(plus::) 

o 

s  .  ass  i  gnc ons ( c  3 , 0 : 1  *  1 ) 

o 

s . f assi gn 

(qlar qlb: : ) 

0 

s . f  ass i gn 

(q2a* q2b: : ) 

o 

s . f mu  1 

(yll/qla*Pa:24,24,24) 

o 

s.  fmul 

(y22,o2a,Pb:24,24,24) 

p 

s . f sub 

(yy,yl 1 , y22 : 24 , 24 , 24 ) 

D 

s • f  sub 

(y»pc»yy:24,24,24) 

D 

s .  f  i  x 

( i  y  f  y :  S ) 

O 

s . anaout 

(iy, 25, -25:8:) 

O 

s .ex i tproc 

(ol us , c  3 : : ) 

O 

s .  var 

(qla:2«) 

D 

s.var 

(q2a:24) 

O 

s .  var 

( i  y :  8 ) 

O 

s . f cons 

(oa,0,0,64, 194:24) 

D 

s  .  f  cons 

(pb,0,0, 192, 193:24) 

O 

s .  f  cons 

(PC/0/0# 1 92 » 128:24) 

O 

s  .end 

(::) 
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B.  OUTPUT  DATA 

1 .  Software  Listing 
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pin  2  =  x61  ;(output) 
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random  access  memory  (lower  half  of  page  255) 
device:  intel  8111-2  1029  bit  (256*9)  static  mos  rami  ic  19 


APPENDIX  C 


ERROR  CORRECTIONS 

The  following  errors  in  the  design  system  were  detected 
and  corrected  during  the  course  of  completing  this  thesis. 

A.  DESIGN  PROGRAM 

The  call  to  subroutine  GETNM  for  subroutine  SYMVAL  was 
corrected.  A  variable  passed  in  the  call,  JJ,  was 
incorrectly  passed  with  a  value  of  zero.  This  was  changed 
to  one  prior  to  the  call  to  GETNM. 

Various  errors  in  branch  statements  existed  throughout 
the  design  program  listing  and  appeared  to  be  the  result  of 
either  an  incorrect  transfer  of  the  original  program  to 
magnetic  tape  or  errors  in  reading  the  magnetic  tape  by  the 
VAX  11/780.  Errors  of  this  nature  were  detected  in  sub¬ 
routines  ARGIF,  SEDUL,  and  TMANAL.  These  have  been 
corrected,  but  validation  of  the  remainder  of  the  program 
must  be  accomplished. 

Subroutine  0UTC0M  writes  the  power  consumption  of  the 
realization  and  its  chip  count  to  the  terminal  screen.  This 
was  observered  to  occur  twice  at  the  end  of  each  execution. 
The  subroutine  was  corrected  so  that  only  one  such  statement 
is  now  displayed  for  the  user.  The  destination  of  the 
duplicate  statement  was  changed  to  the  end  of  the  software 
listing  generated  for  each  successful  realization. 
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Subroutine  SEDUL  calls  a  second  subroutine,  INSERT,  from 
two  possible  locations.  The  first  of  these  used  an  integer 
variable  as  a  parameter  for  a  variable  defined  within  INSERT 
as  real.  This  was  observered  to  prevent  successful 
generation  of  realizations  for  test  input.  The  variable  RZZ 
was  added  as  a  parameter  in  the  calling  routine.  It  is 
equated  to  the  value  of  the  integer  array  element 
0RDER( istop-1 , stop)  prior  to  the  call  to  INSERT. 

INSERT  sets  the  values  of  the  array  ORDER  as  part  of  its 
execution.  The  type  of  ORDER  is  declared  to  be  integer. 
However,  the  values  assigned  to  it  by  INSERT  are  obtained 
from  real  variables.  This  caused  an  integer  overflow  error 
in  the  program  which  resulted  in  termination  of  the  design 
attempt.  The  problem  was  corrected  by  fixing  the  values  of 
the  real  variables  when  equated  to  the  elements  of  array 
ORDER. 

B.  REALIZATION  LIBRARY 

Entry  h.dac  was  missing  the  qualifiers  on  all  skip 
statements.  These  were  determined  and  added. 

Qualifiers  for  skip  statements  in  EPROM  hardware 
primitive  were  set  to  zero.  This  caused  infinite  generation 
of  real  only  memory.  The  qualifiers  were  corrected  to  their 
proper  values. 


Numerous  errors  attributable  to  improper  transfer  of  the 
program  from  magnetic  tape  to  the  VAX  were  detected  and 
corrected.  Validation  of  the  complete  listing  must  be 
performed . 

The  values  of  the  qualifiers  for  the  SKIP  statements  in 
primitives  h.adc  and  h.bufframp  were  incorrect.  These  were 
corrected  to  their  proper  values. 

C.  UNCORRECTED  ERRORS 

Two  errors  are  known  to  exist  in  the  main  program  but 
their  source  could  not  be  determined.  Both  of  these  prevent 
proper  program  execution. 

The  first  error  occurs  when  a  dual  processor  realization 
is  generated.  A  system  error  is  produced  when  the  program 
attempts  to  access  the  file  of  monitor  primitives. 

The  second  error  consists  of  a  failure  to  properly 
evaluate  the  value  of  parameters  enclosed  in  the  symbols  ' < 
>'.  This  prevents  the  proper  utilization  of  library 
primitives  which  use  parameters  as  qualifiers  to  logical 


statements . 
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